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We study possible quantum 17(1) spin liquids in three dimensions with time-reversal symmetry. 
We find a total of 7 families of such 17(1) spin liquids, distinguished by the properties of their emer¬ 
gent electric/magnetic charges. We show how these spin liquids are related to each other. Two of 
these classes admit nontrivial protected surface states which we describe. We show how to access 
all of the 7 spin liquids through slave particle (parton) constructions. We also provide intuitive 
loop gas descriptions of their ground state wave functions. One of these phases is the ‘topolog¬ 
ical Mott insulator’ conventionally described as a topological insulator of an emergent fermionic 
‘spinon’. We show that this phase admits a remarkable dual description as a topological insulator 
of emergent fermionic magnetic monopoles. This results in a new (possibly natural) surface phase 
for the topological Mott insulator and a new slave particle construction. We describe some of the 
continuous quantum phase transitions between the different 17(1) spin liquids. Each of these seven 
families of states admits a finer distinction in terms of their surface properties which we determine 
by combining these spin liquids with symmetry protected topological phases. We discuss lessons for 
materials such as pyrochlore quantum spin ices which may harbor a 17(1) spin liquid. We suggest the 
topological Mott insulator as a possible ground state in some range of parameters for the quantum 
spin ice Hamiltonian. 
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I. INTRODUCTION 

There has been much recent interest in quantum 
spin liquid phases of systems of interacting magnetic 
moments. These phases are fascinating examples of 
ground states characterized by long range quantum en¬ 
tanglement: the corresponding wave functions cannot 
be smoothly deformed into a product state of local de¬ 
grees of freedom. Other examples of such long range 
entangled states include the celebrated fractional quan¬ 
tum Hall states. The structure of the long range entan¬ 
glement dictates the excitation structure of the phase 
in just the same way as the familiar long range order 
(associated with broken symmetry) does in conventional 
ordered phases. In particular a number of unusual ex¬ 
citations - for instance those with fractional quantum 
numbers or statistics, or gapless emergent gauge bosons 
- are possible in such phases. 

This paper is concerned with a particular class of three 
dimensional quantum spin liquid that supports an emer¬ 
gent gapless ‘photon’ as an excitation. It has long been 
recognized [I, 2] that such a photon may be an emergent 
excitation of some underlying physical quantum many 
body system with short range interactions. Specific mi¬ 
croscopic models of quantum phases with emergent pho¬ 
tons were constructed some time ago in Ref. [3-9] in 
diverse systems. In addition to the photon these phases 
support other quasiparticle excitations which couple to 
the photon as ‘electric’ or ‘magnetic’ charge. 

Interest in such phases has been revived following a re¬ 
cent proposed experimental realization[10, II] in certain 
three dimensional pyrochlore oxides. These are materials 
in which there are effective spin-1/2 degrees of freedom 
at the sites of the pyrochlore lattice. A class of such ma¬ 
terials such as Dy 2 Ti 207 or Ho 2 Ti 207 have been stud¬ 
ied extensively both in theory and experiment, and are 
adequately described within the framework of classical 
statistical mechanics [12]. Due to a combination of spin 
anisotropy and exchange interactions, the spins are con¬ 
strained to satisfy an ‘ice rule’ where on each tetrahedron 
of the pyrochlore lattice there are precisely two spins that 
point inward and two that point outward. The low en¬ 
ergy physics takes place within the subspace of states 
satisfying this constraint. Hence these system have been 
dubbed ‘spin ice’. Quantum effects are known to be 
important [13] in a few such pyrochlore magnets which 
have hence been dubbed ‘quantum spin ice’. Examples 
include Yb 2 Ti 207 , Pr 2 Zr 207 , and Tb 2 Ti 207 . In partic¬ 
ular in Yb 2 Ti 207 , the detailed microscopic Hamiltonian 
governing the interaction between the spins has been de¬ 
duced through neutron scattering experiments[10, 14]. 
The parameters of this Hamiltonian are such that quan¬ 
tum effects are surely present and will play a role in 
determining the low temperature physics. 

It is well known[15] that in the spin ice subspace the 


spins form oriented closed loops, and the subspace can 
be parametrized in terms of oriented loop conhgurations. 
Classical spin ice systems are thus described as thermally 
fluctuating loop gases in three dimensions. The loops can 
be viewed as ‘magnetic’ held lines of an artihcial mag¬ 
netic held. Defect conhgurations in the spin ice manifold 
such as a “3-in 1-out” tetrahedron (where 3 spins point 
in instead of 2) correspond to end points of the loops and 
are then identihed with magnetic monopoles[16]. Such 
monopoles have been observed in experiments in the last 
few years[17, 18]. 

In quantum spin ice materials it is natural to expect 
that the physics may be determined by quantum huctua- 
tions of oriented loops. If these loops form a liquid phase 
where the loop line tension is zero the result is a quantum 
spin liquid. This spin liquid supports an emergent gap¬ 
less photon. The associated magnetic held lines are sim¬ 
ply the tensionless magnetic loops. Magnetic monopoles 
(the defect tetrahedra) are now gapped quasiparticle ex¬ 
citations where these held lines end. Ref. [10] proposed 
that this physics may occur in Yb 2 Ti 207 . 

Knowledge of the precise microscopic Hamiltonian for 
Yb 2 Ti 2 07 lends hope for a reliable theoretical assess¬ 
ment of this proposal and for quantitative comparisons 
to experiment. However the microscopic Hamiltonian is 
rather complicated and is hard to solve, either analyt¬ 
ically or numerically. Further as we briehy review (see 
Appendix A) the parameters are such that it is not ob¬ 
vious that it is sufficient to just restrict to the spin ice 
manifold. Finally there is very little global symmetry in 
the model. The only good symmetries are time reversal 
and space group operations. 

What scope is there for theoretical progress in the ab¬ 
sence of reliable methods to study the model Hamilto¬ 
nian? One possibility is to deform the model to a limit 
where it’s ground state may be reliably determined, say, 
by numerical methods. Approximate analytical meth¬ 
ods can then be chosen to reproduce the known result 
in this limit. They can then be extended to the realistic 
model with the hope that they capture the full phase 
diagram. For quantum spin ice such an approach has 
been pursued in Ref. [II] using a slave particle approach 
known as the ‘gauge mean field theory’ (gMFT). A reli¬ 
able limit is provided by considering the XXZ spin-1/2 
pyrochlore model in the Ising limit. This model can be 
studied through quantum Monte Carlo without a sign 
problem and the ground state is a 1/(1) quantum spin 
liquid[8]. Further analytic arguments[5] strongly indi¬ 
cate the structure of the gapped e and m excitations. 
The gMFT correctly reproduces this spin liquid ground 
state. Ref. [11] then extends this slave particle approach 
to obtain an answer for the full phase diagram of the 
model including the parameter regime determined in ex¬ 
periment. This mean held seems to show that the exper¬ 
imentally relevant parameters place the model in a con¬ 
ventional ferromagnetic state rather than a spin liquid. 
However this parameter regime is substantially different 
from the limit where gMFT is known to capture the cor- 
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rect ground state. It is hard to evaluate the accuracy of 
the gMFT prediction for the phase diagram away from 
this limit. In particular other slave particle mean field 
methods are available (for instance Schwinger bosons or 
fermions) and will lead to different phase diagrams. Fur¬ 
ther even when they lead to a 17(1) spin liquid it is not 
clear whether different slave particle methods lead to the 
same phase of matter. 

In this paper, inspired by these developments, we pose 
a different set of questions on which we are able to make 
solid progress. Rather than attempt to solve any partic¬ 
ular microscopic model approximately we constrain the 
general properties of 17(1) quantum spin liquids^ in the 
presence of global symmetries and describe their physics. 
Specifically we focus on time reversal symmetric 17(1) 
quantum spin liquids, and determine the number of dis¬ 
tinct phases and their properties. Time reversal is a 
robust physical symmetry, and the only internal symme¬ 
try in the model describing Yb 2 Ti 207 . We ignore space 
group symmetry both because it simplifies the problem 
and because it is less robust (due to disorder). To further 
simplify the problem we restrict to 17(1) liquid phases 
where the only gapless excitation is the photon. In par¬ 
ticular the magnetic charge (dubbed the M particle) and 
the electric charge (the E particle) are gapped.^ We 
show that there are twenty-two phases which fall into 
seven distinct families of 17(1) spin liquids. The seven 
families of 17(1) spin liquids are distinguished by their 
bulk excitation spectrum, which we tabulate in Table. I. 
Different phases in each family are distinguished by their 
surface states, and one can construct one phase from 
another in the same family by combining with a class 
of phases called symmetry-protected topological states. 
In most parts of this paper except Sec. VIII, we focus 
mainly on the seven families of states, which have clear 
physical differences in the bulk. Therefore we will often 
use the term “phase” instead of “family of phases” when 
the context is clear. Most of the existing microscopic 
models describe only one of these phases which is also 
the one accessed by the gauge mean field theory of Ref. 
[ 11 ]. 

We describe the physics of these seven families of 
states and their interrelationships in many complemen¬ 
tary ways. We show how each of the seven families of 
states may be accessed through slave particle construc¬ 
tions. In some cases we provide more than one slave 
particle construction for the same phase. We describe 
the structure of the distinct ground states in terms of 
distinctions in the loop wave functions. This leads to 
many interesting insights and to predictions for future 
numerical calculations. We determine the properties of 


^ Here we only consider ‘spin liquids’ (or boson liquids) that can 
emerge in the Hilbert space of a purely spin (or boson) system. 

^ We also implicitly assume that, apart from the deconfined 17(1) 
gauge field, there is no other co-existing topological order or 
source of long range entanglement. 


protected surface states that some of these phases have. 
Given these solid results on the possible time reversal 
symmetric 17(1) spin liquids and their properties we may 
ask about how to distinguish them in experiments, and 
about which ones are likely for a particular microscopic 
model. We describe some experimental signatures that 
can help identify which (if any) of these spin liquids is 
realized. We also provide some guides for relating to mi¬ 
croscopic models. A summary of our key results is in 
Sec. II below. 

We emphasize that the distinction between these 
phases is entirely a consequence of the unbroken time 
reversal symmetry. If this symmetry were absent then 
it is possible to go smoothly between any two of these 
phases. The distinction comes from different possible 
implementation of time reversal symmetry. 


Phase 

Electric particle 

Magnetic particle 

El, Mb 

Boson 

Boson 

EbxMb 

Boson, T'^ = — 1 

Boson 

EfMb 

Fermion 

Boson 

EfrMb 

Fermion, = — 1 

Boson 

EbMf 

Boson 

Fermion 

EbrMf 

Boson, T'^ = — 1 

Fermion 

{EfTMf)e 

Fermion, = — 1 

{qm = 2) Fermion 


TABLE I: Families of 17(1) quantum liquids with 
time-reversal symmetry labeled by properties of the 
‘pure’ electric and magnetic charges, qg and qm denote 
electric and magnetic charges, respectively. For the 
electric particle {qe, qm) = (1, 0) and for the magnetic 
particle {qe, qm) = (0,1) except for the last row where 
it is (0,2). For the last phase {EfTMf)g there are more 
fundamental ‘dyonic’ particles which have 
{qe,qm) = (±1/2, ±1), and are bosons. Both the pure 
electric charge and the pure magnetic charge indicated 
in the table can be built up as composites of the dyons 
in this phase. 

Our analysis will be strongly informed by re¬ 
cent progress[19-26] in the theory of interacting gen¬ 
eralizations of three dimensional topological insula¬ 
tors/superconductors (see Ref. [27] for a review of as¬ 
pects directly pertinent to this paper). It is now rec¬ 
ognized that the topological band insulators are special 
examples of a class of quantum states of matter known 
as Symmetry Protected Topological (SPT) phases[28- 
30]. These states are seemingly conventional in the bulk 
- they are gapped and have no exotic excitations but 
nevertheless have non-trivial surface states that are pro¬ 
tected by symmetry. But what role do they play in de¬ 
scribing the quantum spin liquids of interest in this pa¬ 
per? The answer is that starting with one kind of 17(1) 
spin liquid we may generate others by putting one of the 
emergent quasiparticles {E or M) into an SPT state. 
For 17(1) quantum spin liquids this point of view was 
initiated in a previous paper[20] by the present authors. 
Ref. [20] considered SPT states of bosonic particles, and 
demonstrated their utility in understanding some 17(1) 
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quantum spin liquids. This point of view will be fully de¬ 
veloped in the present paper and lead to a complete and 
more insightful description of all time reversal invariant 
C/(l) spin liquids with gapped matter. In particular we 
will exploit recent exciting developments on fermionic 
SPT states [23-26] that were not understood when Ref. 
[20] was published to obtain this complete picture. 

It is important to point out that there is no one-to- 
one mapping between SPT phases with global U{1) and 
time reversal symmetries, and 17(1) quantum spin liquids 
with time reversal. They both have different classihca- 
tions. For example, we will show that two different SPT 
states reduce after gauging to the same physical 17(1) 
spin liquid. 

II. SUMMARY OF RESULTS 

Here we briefly summarize some of our key results. 
This section will also serve as an outline for the rest of 
the paper. 

1. We first establish that there are 7 distinct fami¬ 
lies of time reversal invariant 17(1) liquid phases 
in 3(7 distinguished by their bulk spectra in Sec¬ 
tions III, IV and V. A partial characterization of 
these phases is obtained by plotting the spectrum 
of emergent quasiparticles - the charge-monopole 
lattice - in the 17(1) gauge theory. We show that 
six of these 7 families have the charge-monopole 
lattice of Fig. 1 while the remaining one has the 
charge-monopole lattice of Fig. 2. We provide a 
first cut description of these 7 families of phases in 
these sections and relate them to existing construc¬ 
tions of 17(1) liquids. One of these phases (dubbed 
EbMb) is the state accessed by gMFT while some 
others are states accessed by Schwinger boson or 
Abrikosov fermion constructions. The unique fam¬ 
ily described by Fig. 2 includes the so-called ‘Topo¬ 
logical Mott Insulator’ discussed in Ref. [31]. For 
reasons described later this is denoted (Fl/^My)^ 
in this paper. The family denoted E^TMf has not 
been described explicitly in the literature. 

2. We describe how these phases are related to each 
other in Sec. VI. This is enabled by recent 
advances in our understanding of SPT phases of 
bosons/fermions with global 17(1) and time rever¬ 
sal symmetries. We continue the point of view 
adopted in our previous work[20] showing that 
given one 17(1) liquid we can obtain others from 
it by putting one of the emergent quasiparticles in 
an SPT phase. We are thus able to obtain a rather 
complete understanding of how these 7 phases are 
related to each other. 

3. The conventional description of the {EfTMf)g (the 
‘Topological Mott Insulator’) is that it is a topo¬ 
logical insulator formed by emergent fermionic 


Kramers doublet spinons that are coupled to the 
17(1) gauge field as electric charge. We show in 
Sec. VI that this phase has a remarkable dual 
description as a topological insulator of emergent 
fermionic magnetic monopoles. 

4. We discuss the possibility of protected surface 
states at the interface with the vacuum for these 
spin liquids in Sec. VII. In subsection VIIA we de¬ 
scribe criteria that determine when such protected 
surface states will form. We argue that precisely 
two of the 7 families {{EfTMf)g and Ei^-^Mf) are 
required to have a non-trivial surface state. The 
possibility of a surface spinon Dirac cone for the 
{EfTMf)g is well known[31]. The dual description 
of this phase as a monopole topological insulator 
naturally leads to an alternate possible ‘dual’ sur¬ 
face state where there are an odd number of gap¬ 
less monopole Dirac cones (and no spinon Dirac 
cone). We then describe the surface of E^TMf - 
the simplest possibility is to have two monopole 
Dirac cones. 

5. For any given bulk spectrum there can be more 
than one phase corresponding to distinct surface 
properties. When these are taken into account we 
find a total of 22 distinct phases. These are ob¬ 
tained from the 7 basic phases by combining them 
with SPT phases of bosons/spins protected by time 
reversal alone. Interestingly, in some cases, the 
spin liquid can ‘absorb’ an SPT phase so that the 
combination is not in a distinct phase. In other 
words not all T-reversal symmetric SPT phase stay 
distinct from trivial phases when combined with a 
spin liquid. Similar phenomenon also appears in 
two dimensional Z 2 spin liquids [32]. 

6. We show, in Sec. IX how to access all the 7 ba¬ 
sic phases through parton constructions on spin 
models. In particular we show how the standard 
fermionic parton construction of spin-1/2 systems 
enables access to 5 of the 7 phases (the excep¬ 
tions being E^Mb described by gMFT and the 
EbrMb described by Schwinger bosons). For the 
topological Mott insulator {EfTMf)g we describe 
a dual parton construction in terms of monopoles 
that is distinct from the conventional one in terms 
of spinons. With a view toward obtaining input 
on microscopies, we obtain some no-go results on 
these parton constructions in Sec. IX A if the phys¬ 
ical system consists of Kramers doublet spin-1/2 
degrees of freedom. 

7. We provide an intuitive physical picture of the 
ground state wave function for these spin liquids 
in terms of fluctuating loop configurations of elec¬ 
tric or magnetic field lines in Sec. X. 

8. We describe some of the remarkable continuous 
quantum phase transitions between these differ- 
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ent spin liquids. We particularly focus on phase 
transitions out of the topological Mott insulator 
{EfTMf)g. We provide a theory for a second or¬ 
der transition from this phase to others where the 
electric charge is a boson (either Kramers singlet 
or doublet). We also provide a theory for a dif¬ 
ferent second order phase transition between two 
different phases where the electric charge changes 
from Kramers doublet to Kramers singlet. 

9. In Sec. XIV we consider the relevance of these re¬ 
sults to current and future realizations of C/(I) spin 
liquids in experimental systems. For pyrochlore 
spin ices based on Kramers doublet spin systems 
we discuss the possible 17(1) spin liquids that may 
obtain. Apart from the one suggested by gMFT, 
we argue, based on the parton construction, that 
the other natural candidate is the topological Mott 
insulator (EfTMf)g. A strong coupling expansion 
of the lattice parton Hamiltonian coupled to the 
U{1) gauge field yields at leading order a spin 
Hamiltonian of the form appropriate to the py¬ 
rochlore spin ices but with parameters different 
from the ones where gMFT is expected to be re¬ 
liable. We suggest on this basis that some py¬ 
rochlore spin ices may be in the topological Mott 
insulator phase. 

We also describe some experimental distinctions 
between these phases which may be useful in iden¬ 
tifying them. 

Several appendices contain peripheral details. 

III. PRELIMINARIES 

We begin with some simple but powerful observations. 
We are interested in time reversal symmetric 17(1) liq¬ 
uids of spins/bosons in which the only gapless excitation 
is the photon. To distinguish different phases it is ap¬ 
propriate to focus on the gapped emergent quasiparticles 
that couple as electric or magnetic charges to the pho¬ 
ton. Time reversal symmetry constrains the possibilities 
in many important ways as we now describe. 

A. Charge-monopole lattice 

We denote the electric charge qe and magnetic charge 
qm- We use notation in which the total electric flux 
is Airqe and the total magnetic flux is 2'Kqra- To be 
general we must allow for the most fundamental emer¬ 
gent particles to be “dyons”, i.e particles that carry 
both electric charge and magnetic charge. For any 
pair of dyons with charges {qe = Qe,qm = Qm) and 
(<Ze = QejQm = Q'm) there is a generalized Dirac quanti¬ 
zation condition[33, 34]: 

QeQ'm. - QniQ'e = (1) 


where n is an integer. 

For each particle with charges {Qe,Qm) there will be 
an antiparticle with charges {—Qe, —Qm)- Note that the 
particle and antiparticle automatically satisfy the Dirac 
quantization condition. We will use the natural conven¬ 
tion that under time reversal the magnetic fields are odd 
and the electric fields are even. Then for any particle 
with charges {Qe,Qm) there is a time reversed partner 
with charges {Qe, —Qm)- Applying the Dirac quantiza¬ 
tion condition to these two particles we obtain the re¬ 
striction 

2QeQm = integer ( 2 ) 

By combining {Qe,Qm) with {Qe, —Qm) we can pro¬ 
duce a particle that is a pure electric charge {2Qe,0)- 
Similarly by combining {Qe,Qm) with {—Qe,Qm) (the 
antiparticle of the time reversed partner) we obtain a 
pure magnetic charge (0, 2Qm). Thus time reversal in¬ 
variance guarantees that there are always both pure elec¬ 
tric and pure magnetic charges in the theory. 

Consider the smallest pure electric charge. We choose 
units in which this has qe = 1 (and by definition has 
qm = 0). Let the smallest pure magnetic charge have 
strength g (and qe = 0). Applying the Dirac condition 
to the pure electric charge and the pure magnetic charge 
we get 

g = integer (3) 

If there are no other restrictions the smallest allowed g 
is I. As is well known the Dirac condition requires that 
pure electric and charges are quantized to be integers 
(in our units). If there are dyons with charges {Qe,Qm) 
it follows that 2Qe must be an integer. Thus we have 
two basic possibilities Qe = 1 or Qe = 5 ■ In the former 
case there are no further restrictions on g beyond Eq. (3) 
and we have g = 1- In the latter case we can apply Dirac 
quantization to the ( 0 , g) and {\,Qm) particles to obtain 

g = 2 X integer (4) 

Thus if there are charge-I/2 dyons the minimum pure 
magnetic charge is 2. Further we must have Qm = 1 for 
the charge-I/2 dyon. 

We thus have two classes of possible states which 
are distinguished by the geometry of the lattice of al¬ 
lowed charges and monopoles. In one class the charge- 
monopole lattice is as shown in Fig. I. Here all emergent 
quasiparticle excitations are obtained from two elemen¬ 
tary quasiparticles - the E particle with {qe,qm) = (1)0) 

and the M particle with {qe,qm) = (0,1). In the sec¬ 

ond class the charge-monopole lattice is shown in Fig. 
2. Here the full set of emergent particles can be built 
out of the two dyons with {qe,qm) = (i|, I). 

These two general possibilities can also formally be 
distinguished in terms of the low energy effective La- 
grangian for the photon after integrating out the E and 
M particles. This takes the form 

^eff — ^Max T ^9 (b) 
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FIG. 1: Charge-monopole lattice at 0 = titt with n even. 



FIG. 2: Gharge-monopole lattice at 6 * = titt with n odd. 


The first term is the usual Maxwell term and the second 
is the ‘theta’ term: 

= ( 6 ) 

where E and B are the electric and magnetic fields re¬ 
spectively. 

As is well known time reversal restricts the allowed 
values to 9 = nn with n an integer, n even corresponds 
to Fig. I and n odd to Fig. 2[35]. 

Each of these two charge-monopole lattices can po¬ 
tentially be realized in several ways depending on the 
statistics of the quasiparticles and transformation under 
time reversal. We next describe the constraints on these. 


B. Quasiparticle statistics and symmetry 
realization 

For a spin or boson system, in the microscopic Hilbert 
space, excitations created by local physical operators 
must clearly be bosonic. The emergent quasiparticles 
E or M are however not created by local operators. For 
instance to create an E particle it is also necessary to 
create the ‘electric’ field lines that emanate from it and 
which extend out to arbitrary long distance. A formal 
way of describing this is to say that the ‘creation’ oper¬ 
ator for E (or M) alone is not gauge invariant. Creating 
E without the associated electric field violates Gauss law 
and hence is not in the physical Hilbert space. As E and 
M are not created by local physical operators there is 


no restriction that they must be bosonic. In three space 
dimensions they can thus be either bosons or fermions. 
Very recently it has been shown[24, 36] however that in a 
strictly 3d spin/bosonic system (as opposed to systems 
that can only appear the boundary of a 4 -|- 1 dimen¬ 
sional system) that E and M cannot simultaneously be 
fermionic. We will return to this below. 

Time reversal symmetry acts in a simple way on phys¬ 
ical states in the Hilbert space of spins/bosons. The 
time reversal operator (T) is anti-unitary and satisfies 
= -1-1 on all physical states. This should be con¬ 
trasted with electronic systems where = — 1 for an 
odd number of electrons which leads to Kramers degen¬ 
eracy. Let us now discuss the possible action of time 
reversal on the emergent E and M particles. Quite gen¬ 
erally the structure of the emergent Maxwell equations 
implies that the electric charge is even while the mag¬ 
netic charge is odd under time reversal^. Thus the E 
particle and its time reversed partner TE only differ by 
a local operator. Then acting on the E particle has a 
well defined value. Now as E itself is not local, it could 
have = — 1 and hence be a Kramers doublet (we re¬ 
view some more details in Appendix B). In contrast M 
and TM do not differ by a local operator. Then there 
is no meaning to asking whether M is Kramers or not. 
Specifically acting on M can be shifted by a gauge 
transformation to have any value[20]. 

Finally though E or M may be a fermion, and E may 
be a Kramers doublet, composite excitations formed out 
of them that carry zero electric and magnetic charge are 
physical excitations, and hence must be bosonic Kramers 
singlets. 

Starting with these simple but powerful observations 
we proceed to describe all the distinct time reversal in¬ 
variant C/(l) spin liquids where the photon is the only 
gapless excitation. 

IV. PHASES WITH 6 = 0 

We first describe phases in which the parameter 9 = 0, 
i.e the charge-monopole lattice is given by Fig. 1. Here 
we distinguish two broad classes of phases depending 
on whether the M particle is a boson or fermion. We 
describe each in turn. 

A. Bosonic monopole 

It is clear first that there are 4 distinct phases where 
M is a boson. The E particle may either be a bo- 


® Strictly speaking what we call electric and what we call magnetic 
is a matter of convention; the 1/(1) gauge theory is self-dual so 
that we can interchange the definition of E and M. Maxwell’s 
equations imply that the electric and magnetic charge trans¬ 
form oppositely under time reversal. It is natural to adopt the 
convention that the magnetic charge is time reversal odd. 
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son or a fermion and be Kramers singlet or doublet. 
Let us understand better these 4 phases. We label 
them EfMh, EfTMf,, respectively, with 

the subscripts b, f describing the statistics, and the sym¬ 
bol T referring to Kramers degeneracy. Some of these are 
obtained through familiar constructions. 

The EbMb phase is the one constructed in most of the 
existing microscopic models[3, 5-9]. It is also the state 
accessed by the gauge mean field theory of Ref. [11]. If 
in addition to time reversal there is a global 1/(1) sym¬ 
metry, an intuitive way to understand this phase was 
described in Ref. [7] by obtaining it from a proximate 
long range ordered phase {i.e with broken 1/(1) symme¬ 
try) through proliferating appropriate vortex loops. 

The phases EhxMh, EfxM^ are accessed by the stan¬ 
dard Schwinger boson or Abrikosov fermion representa¬ 
tion of the physical spin. It is well known that in Sd 
these representations can lead to stable 1/(1) spin liquid 
phases with gapped electric and magnetic charges. 

For spin systems with spin rotation symmetry, it is 
instructive to obtain Ei,xMi, by starting with a semi- 
classical description of a Neel antiferromagnet as follows. 
Consider a collinear Neel state of a quantum antiferro¬ 
magnet in 3d. The corresponding order parameter man¬ 
ifold is An effective field theory description of the 
long wavelength fluctuations of the Neel order parameter 
is provided by the quantum non-linear sigma model in 
3-1-1 space-time dimensions with the Euclidean action: 


SnLuM 


1 


d^xdr 




(7) 


Here n is the local orientation of the Neel vector, and 
c is the spin wave velocity. In 3d the order parame¬ 
ter manifold allows for point defects known as ‘hedge¬ 
hogs’ corresponding to n 2 (S'^) = Z. This Neel state 
may be quantum disordered without proliferating these 
hedgehogs. A convenient framework to describe this is 
through a CP^ representation: n = z^'^az where z is 
a two-component complex spinor. Importantly under 
time reversal n —> —n and Za —t Thus z is 

a Kramers doublet. The z representation introduces a 
1/(1) gauge redundancy {z{x,t) —>■ z{x,t)). The 

sigma model action represented in terms of z naturally 
includes a compact C/(l) gauge field a^. It is well known 
that the monopoles of the correspond, in the Neel 
ordered state, to the hedgehogs of the n field. Quantum 
disordering the Neel state corresponds to gapping out the 
z-particles. If in addition the monopoles stay gapped the 
result is precisely a 1/(1) spin liquid. Further the z get 
identihed with the E particle and the M with the rem¬ 
nants of the hedgehog. Clearly A is a Kramers boson. 
In the semiclassical limit the hedgehog is also a boson 
and consequently so is the M particle in the 1/(1) spin 
liquid. Thus the phase we obtain is precisely the E^rMb 
1 /(1) spin liquid. 

Finally, a microscopic model for the EfMi, phase was 
constructed in Ref. [37]. In Sec. VI we describe how it is 


related to the other phases, in particular to the simple 
EbMf, phase. 


B. Fermionic monopole 

We now consider cases in which the M particle carries 
fermion statistics. If M is a fermion, recent work[24, 36] 
shows that the E particle cannot also be a fermion in 
a strictly three dimensional system. With a bosonic E 
particle there are however still two distinct possibilities 
corresponding to whether it has = -1-1 or = — 1, 
i.e whether it is a Kramers singlet or doublet. In obvious 
notation we label these two phases E^Mf and E^rMf. 

We show in Section IX how to access these phases 
through a parton construction. 


V. PHASES WITH 9 = tv: “TOPOLOGICAL 
MOTT INSULATOR” 

We now discuss time reversal symmetric 1/(1) spin liq¬ 
uids with 9 = TT. We will see that there is precisely one 
such phase. 

First let us discuss the statistics of the elementary 
dyons with charges {qe,qm) = (^jTl)- We note that 
these are interchanged under time reversal. Thus they 
are required to have the same statistics, be they are both 
bosons or both fermions. However we can argue that 
they cannot both be fermions. To see this most simply 
we note that the ((|, 1) and the (|, — 1) dyon are relative 
monopoles, i.e each one sees the other the way an electric 
charge sees a monopole. If they were both fermionic we 
would have a realization of the “all-fermion” 1/(1) gauge 
theory in a strictly 3-1-1 dimensional system which we 
know is not possible [24, 36]. Therefore we conclude that 
both these dyons must be bosons. 

Now consider the bound state of these two dyons. As 
this has qe = l,qm = 0 we identify it with the ‘elemen¬ 
tary’ pure electric charge in this phase. Precisely this 
bound state was analysed recently in Refs. [21, 24] while 
studying correlated topological insulators, and shown 
to be a fermion with = —1, i.e, a Kramers dou¬ 
blet. In brief these two dyons see each other as relative 
monopoles. This leads to the Fermi statistics of their 
bound state. The Kramers degeneracy can be simply 
understood by first calculating the angular momentum 
of the 1/(1) gauge field. It is readily seen that this is 
quantized to be 1/2. Combining this with the observa¬ 
tion that time reversal inverts the relative coordinate of 
the two dyons leads to = — 1 for their bound state. 
Thus the statistics and symmetry properties of the ele¬ 
mentary electric charge are uniquely determined for this 
charge-monopole lattice. 

Next consider the elementary pure magnetic charge 
which has qe = 0,qm = 2. This can be obtained as the 
bound state of the ((|,1) and ((—/,1) dyons. These 
are also relative monopoles and hence their bound state 



is a fermion. Now time reversal does not interchange 
these two dyons and hence the argument above for the 
Kramers structure of the pure electric charge does not 
apply. This is of course in line with the earlier argument 
that it is meaningless to ask if M particles are Kramers 
or not. 

We thus see that the structure of both the elemen¬ 
tary electric charge and the elementary magnetic charge 
are uniquely determined for this charge-monopole lat¬ 
tice. In addition the statistics and symmetry properties 
of the elementary dyons is also fixed. Thus there is pre¬ 
cisely one time reversal symmetric 17(1) spin liquid phase 
corresponding to 9 = tt. Given these properties of the 
elementary pure electric and magnetic charges we denote 
this phase The subscript 0 is a reminder that 

these pure charges are composites of more fundamental 
dyons. 

This phase may be constructed within slave particle 
methods. Let us begin with EfTMi, where the E particle 
is a Kramers doublet fermion, and has a conserved elec¬ 
tric charge that is even under time reversal. We then put 
this E particle into a topological band insulator phase. 
It is well known[38] that the topological band structure 
leads to a 0 = TT term in the action for a 17(1) gauge 
field that couples to the E particle. This slave particle 
construction of the {EfTMf)g 17(1) spin liquid was dis¬ 
cussed in Ref. [31] and dubbed the ‘topological Mott 
insulator’. Ref. [31] also suggested that this phase may 
be realized in Y 2 lr 20 j though this has turned out to be 
unlikely. 

Later we will describe a completely different slave par¬ 
ticle construction of this phase. 


VI. RELATIONSHIP BETWEEN THE PHASES 

We have thus completed the description of Table 1. 
We now describe how these seven different phases are re¬ 
lated to each other. In our previous work (Ref. [20]) we 
addressed this for a subset of these phases, and showed 
that they can be related to different Symmetry Protected 
Topological (SPT) phases[30] of one of the emergent ex¬ 
citations, similar to what has been discussed for topolog¬ 
ical orders [39]. Here we will continue this point of view 
to develop a detailed understanding of the relationship 
between all seven phases, which is summarized in Fig. 3. 
This exercise adds much new insight, and provides for 
new constructions of some of these phases. It also helps 
us obtain theories for some of the quantum phase tran¬ 
sitions between these spin liquids. 

We have already discussed how {EfTMf)g may be un¬ 
derstood as a topological insulator of the EfT particle. 
So we now turn to the other phases. 

Let us start from the simple E^Mb phase which can be 
obtained straightforwardly in microscopic models. We 
can obtain new phases by either putting M or i? in 
bosonic topological insulator phases. E transforms under 
Ueg{l) x Z 2 (meaning that the electric charge is T-even) 



FIG. 3: Relationship between different 17(1) spin 
liquids. Two phases connected through a line share a 
common fundamental particle {E or M), and can be 
viewed as different SPT phases formed by the common 
particle. In Sec. XI we describe some intersting 
continuous phase transitions between the phases 
connected through thick red lines. 


while M transforms under Umg x Z 2 (meaning that the 
magnetic charge is T-odd) where 17eg(l) is the electric 
gauge transformation and Umg{^) is the magnetic gauge 
transformation. 

We discuss this first for the M particle. Con¬ 
sider bosonic topological insulators with global symmetry 
C/(l) X Z 2 ■ There are a total of 16 such phases corre¬ 
sponding to classification by the group Z^. These can be 
obtained from 4 ‘root’ phases (the 4 generators of Z|) 
and taking their combinations. Two of these root phases 
are protected by time reversal alone while the remaining 
two require the full 1/(1) x Z 2 symmetry. Now consider 
coupling these bosons to a dynamical 1/(1) gauge field, 
i.e gauging the global 17(1) symmetry. The two root 
phases whose distinction requires also the C/(l) subgroup 
then potentially lead to gauge theories with distinct bulk 
excitations^. Taken together with their combinations we 
get a total of four potentially distinct 1/(1) spin liquids. 
The understanding of such bosonic SPT phases shows 
that these are precisely the four U{1) spin liquids with a 
bosonic monopole {EbMb, EfMb, EbrMb, E/rMb) dis¬ 
cussed in Sec. IV A. 

Next consider starting with EbMb and putting E in 
a boson topological insulator. Such insulators with 
1/(1) X Z2 symmetry are classified by with 3 root 
phases. Of these only one is protected by the full 
1/(1) X Z2 symmetry. Coupling the E, when it forms this 
SPT state, to a dynamical 1/(1) gauge field then leads to 
fermion statistics of the M particle[20, 21]. Thus we ob¬ 
tain the EbMf spin liquid. The fermionic statistics of the 
M particle can be understood from a 0-term in the gauge 
theory with 0 = 27r (Ref. [19]) and was called the “statis¬ 
tical Witten effect” [21]. Likewise starting with EbrMb 
state one can also put the Kramers bosonic charge E in 
an SPT state and obtain the EbxMf state. 


^ The other phases correspond to combining the P(l) liquids with 
SPT paramagnets protected by time reversal alone. We defer a 
discussion of these to Sec. VIII 
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Let us now understand the phases obtained by start¬ 
ing with EfjMf and putting M in a topological insulating 
phase. We first recall that the monopole transforms un¬ 
der 17mg(l) X Z 2 where the Umg{^) is the (magnetic) 
gauge transformation. Free fermions with global sym¬ 
metry 17(1) X Z 2 can form topological band structure 
classified by Z, i.e there are distinct phases indexed by 
an integer n which counts the number of Dirac cones 
at the surface. With interactions this collapses to a Zs 
classification[25, 26]. Hence we only need to consider 
n{mod8). Of these n = 4 is protected by Z 2 alone. We 
now argue that if the global 17(1) is gauged, as appro¬ 
priate in the 17(1) spin liquid, then n = 0,2 (and only 
these) lead to distinct (at the level of bulk excitations) 
phases 

n = 0 corresponds simply to the E^Mf phase. Inter¬ 
estingly n = 2 corresponds to the E^rMf phase[25]. 


A. A puzzle 

When n = 1, the 1/(1) gauge field acquires a 9 term 
at 0 = TT. This is an example of a ‘topological Mott 
insulator’ that seems distinct from the one discussed in 
the previous section. In contrast to the description of the 
{EfTMf)g as a topological insulator of EfT, here the 9 
term originates from the M sector and leads to a ‘dual’ 
Witten effect whereby the E particle acquires magnetic 
charge 1/2. How do we reconcile this with our claim that 
the list of 7 phases is complete? 


B. Resolution: a dual description of the 
Topological Mott Insulator 

The resolution of the puzzle above is that the n = 1 
topological insulator formed by the Mf particles is actu¬ 
ally identical to the {EfTMf)g phase. To see this con¬ 
sider the charge-monopole lattice of the n = 1 monopole 
topological insulator. This is shown in Fig. 4. Clearly it 
is is very similar to that of {EfTMf)g. In Sec. HIA we 
made the choice of units that the minimum pure elec¬ 
tric charge is 1. This leads - when there were funda¬ 
mental dyons - to a minimum pure magnetic charge of 
2. We could equally well have chosen units so that the 
minimum pure magnetic charge is 1. Dirac quantization 
(together with the existence of fundamental dyons) then 
would demand that the minimum pure electric charge 
is 2. This corresponds to taking the lattice in Fig. 2, 
shrinking the g^-axis by a factor of 2, while expanding 
the ge-axis by a factor of 2. This converts the lattice 
in Fig. 2 to that in Fig. 4. Clearly this change in the 
unit choice does not change the physics. In particular 
we correctly obtain that the pure electric charge (which 
has charge-2 in these units) is a Kramers fermion. 

Thus this is not a new phase but rather is included in 
our list. 



FIG. 4: Charge-monopole lattice obtained by gauging 
the n = 1 Mf topological insulator. It is identical to 
Fig. 2 after rescaling the two axes as explained in the 
text. 


Let us revert back to the units where the minimum 
pure electric charge is 1. We see that remarkably there 
are two equivalent descriptions of the {EfTMf)g phase 
which are dual to each other. We can either describe it as 
a topological insulator of A/t or as a n = 1 topological 
insulator of My. This leads to a number of interesting 
consequences which will be explored in subsequent sec¬ 
tions. 

In Appendix C we show why My topological insula¬ 
tors with other values of n do not lead to distinct U{1) 
spin liquids. We note, and will discuss in greater detail 
below, that a duality similar to the one above also exists 
for the spin liquid A^rMy. First it can be thought of as 
the n = 2 U{1) xT topological insulator of the Mf par¬ 
ticle. Equivalently it can also be viewed as a boson topo¬ 
logical insulator with Kramers charge (the Ei,t)- This 
duality will have interesting consequence for the surface 
state, which we discuss in Sec. VH B, and for a loop wave 
function for this phase (Sec. XA. 


VII. SURFACE STATES 

The understanding of the connection between these 
17(1) liquids and SPT states immediately raises the ques¬ 
tion of whether there are non-trivial surface states at 
the boundary between any of these spin liquids and 
the vacuum. In this section we discuss the necessity 
(or lack there of) of nontrivial surface states on general 
grounds. We then discuss two interesting examples: the 
{EfTMf)g and the E^xMf spin liquids, which necessar¬ 
ily have non-trivial surface states. In Appendix D we 
discuss surface states of the other phases, and the inter¬ 
face between different phases. 


A. Why surface states? 

To see why (or why not) there should be a surface state 
between a 1/(1) spin liquid and the vacuum, we should 
first understand what exactly is a vacuum. Since the 
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Trivial 

= = • 

Vacuum 


FIG. 5: The “wall” between a 17(1) spin liquid and a 
Higgsed vacuum. A particle {E or M) tunnels through 
the wall and becomes a trivial boson, which 
subsequently condenses and forms the vacuum. Some 
nontrivial excitation must be left behind on the wall 
after the tunneling process. 


17(1) gauge field “disappears” in the vacuum, we should 
really think of the vacuum as a confined phase of the 
U{1) gauge theory. In (3 + 1) dimensions, a 17(1) gauge 
theory can confine in two ways: by condensing (Higgs- 
ing) either the E or M particle. Therefore if either the E 
or M particle is a non-Kramers boson and the vacuum is 
simply the condensate of that particle, the surface state 
will be featureless. However, if either the E or M par¬ 
ticle carries nontrivial quantum number (fermion statis¬ 
tics or Kramers degeneracy), it cannot directly condense 
and form the vacuum. Instead, it should go through a 
“wall” that converts it into a trivial boson, so that the 
trivial boson could condense and form the vacuum. The 
“wall” then forms the surface between the 17(1) spin liq¬ 
uid and the vacuum, and obviously something nontrivial 
is needed on the wall for the conversion (Fig. 5). A sim¬ 
ilar reasoning was used in a slightly different context in 
Ref. [21]. We name the A-particle-converting wall as 
A-wall, and likewise M-wall for the M-converting wall 
(a similar notation was also used for Z 2 spin liquids in 
2D, for example in Ref. [40]). Since the A-condensate 
and the M-condensate are really the same vacuum, the 
A-wall and M-wall can evolve into each other through 
phase transitions on the surface, without actually chang¬ 
ing the vacuum. 


B. Surface of (A/tM/)^ 


As a particularly interesting example, we now de¬ 
scribe the surface of the 9 = tt spin liquid. Within 
the Abrikosov fermion slave particle construction this 
state (see Ref. [31] ) appears as a topological insulator 
of the EfT particle. A natural conclusion is that the 
surface will have an odd number of Dirac cones of the 
EfT particle which are then coupled to the bulk gapless 
17(1) gauge field. In the language of Sec. VII A, this is 
an M-wall, since a charge-1/2 dyon can tunnel through 
the wall and become a pure monopole M, which can 
subsequently condense and form the vacuum. The M- 
converting phenomenon is simply a manifestation of the 
parity anomaly[41] of the surface Dirac cone: a tunneling 
process that changes the surface flux by 27r will change 
the total charge[42] on the surface by ±1/2. 

The surface could, just like for the topological insula¬ 
tor, break time reversal symmetry. The Dirac fermion 
will then be gapped, but the surface will have a half¬ 
integer hall conductance for the gauge field. Alternately 
the surface could spontaneously enter a Higgs phase by 
condensing a bosonic Cooper pair formed out of the Ejt 
particles. This corresponds to the surface superconduc¬ 
tor, in the underlying topological insulator. Again the 
matter fields are gapped on the surface, but there will 
be vortex excitations with nontrivial fusion and braiding 
properties. Finally for the underlying topological insu¬ 
lator gapped symmetry preserving surface states with 
anyonic excitations are also possible[43-46]. In the 17(1) 
spin liquid these anyons will, if charged, be coupled to 
the 17(1) gauge field. In all these cases, the surface states 
are M-walls. 

The alternate view of this state as a monopole topolog¬ 
ical insulator leads naturally to a very different gapless 
symmetry preserving surface state: The n = 1 Mf topo¬ 
logical insulator has at the surface a single Dirac cone 
formed out of the Mf particles. This Dirac cone is neces¬ 
sarily at the neutrality point as the density oi Mf is odd 
under time reversal. In contrast to the electric Dirac 
cone, tunneling a dyon through the wall gives a pure 
electric charge. Therefore this monopole Dirac cone is 
an A-wall, and is a very different ‘dual’ possibility for 
the surface state of the topological Mott insulator. 


Therefore for the two phases {EfTMf)g and E^xMf, 
in which all the fundamental particles are nontrivial, a 
nontrivial surface state is necessary, no matter how we 
view the vacuum. For the other five phases in Table I, at 
least one of the E and M particle is a trivial boson, hence 
the surface is allowed to be featureless. However, they 
can nevertheless have nontrivial surface states if we view 
the vacuum differently. This will be particularly relevant 
if we consider the interface between two different C/(l) 
spin liquids, which we discuss in Appendix D 


Given a realization of the {EfTMf)g phase which of 
these many surface phases is realized will be determined 
by microscopic details. As the parameters of a micro¬ 
scopic Hamiltonian are tuned while keeping the bulk in 
this phase the surface may undergo phase transitions be¬ 
tween these various phases. In particular the Efx Dirac 
cone state may transition into the M f Dirac cone as the 
parameters are varied. This is related to a “dual” Dirac 
liquid state that can be realized on the surface of a topo¬ 
logical insulator[47, 48]. 
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C. Surface of EtrMf 

The EbTMf C/(l) spin liquid is obtained naturally as 
a n = 2 Mf topological insulator. This point of view 
immediately tells us that the surface will have a gapless 
state with 2 Dirac cones (at neutrality) of the Mf par¬ 
ticles. These will be coupled to the bulk 17(1) gauge 
field. In this case all Mf topological insulators with 
n = 2 {modA) lead to the same EhrMf 17(1) phase in the 
presence of the bulk gauge field. Therefore the gapless 
symmetric surface state may have n = 2 (mod 4) num¬ 
ber of Dirac cones. This surface is an if-wall since it 
converts the Kramers E particle into a non-Kramers one 
upon tunneling. Just as in the previous subsection alter¬ 
nate surface states with deconhned anyonic excitations 
or breaking symmetries are also possible. 

The EbxMf state can also have a nontrivial M-wall. 
We describe it in more detail here as an example of a sur¬ 
face state with gapped matter fields. It will also be use¬ 
ful in constructing a loop wavefunction of the spin liquid 
phase, which we discuss in Sec. X A. In such a surface the 
electric charge is gapped and form a Z 2 topological order, 
with topological quasi-particles denoted as {l,e, m, e}. 
The symmetries are assigned to the quasi-particles as fol¬ 
lows: e and m carry electric charge qe = 1/2, T switches 
e and m, and e is charge-neutral and T-singlet. Because 
e and m have a mutual 7r-statistics and are exchanged 
under E, the bound state em will have = — 1. There¬ 
fore = e{mt) also has = —1, and likewise for rn?, 
which is consistent with the bulk physics since they cor¬ 
responds to the charge-1 boson. Following the logic of 
Ref. [21], tunneling the M through this wall will change 
its statistics from fermion to boson. 

This surface state is precisely the surface of a bo¬ 
son topological insulator formed by the Ei,t particle. 
Such topological insulators of Kramers doublet bosons 
have not been discussed much in the literature (as far as 
we know). However their physics is easily deduced us¬ 
ing methods developed[19-21] for non-Kramers bosonic 
topological insulators. Apart from the surface topologi¬ 
cal order described in the previous paragraph, just like 
their non-Kramers cousins, the E^t topological insula¬ 
tor haa a 0 = 2 tt response, and associated fermionic 
monopoles. 


VIII. COMBINING 1/(1) SPIN LIQUIDS WITH 
TOPOLOGICAL PARAMAGNETS 

So far we have identified two phases with the same 
bulk excitation spectrum as the same phase. This is rea¬ 
sonable if the relevant experimental probes are only de¬ 
tecting the bulk physics. However, additional structure 
needs to be considered if one is also interested in the sur¬ 
face states of the 17(1) spin liquids. In particular, one can 
combine a 17(1) spin liquid with a symmetry-protected 


topological (SPT) state^. This does not change the bulk 
spectrum in any nontrivial way, but may produce dis¬ 
tinct surface states. 

SPT states with only time-reversal symmetry are also 
called “topological paramagnets” [19]. In three dimen¬ 
sions, it is known that there are three nontrivial topolog¬ 
ical paramagnets. Together with the trivial state, they 
form a structure, which simply means combining two 
copies of the same state always produces a trivial state, 
and combining two distinct states gives the third distinct 
state. In the notation of Ref. [20], the three nontrivial 
states are labeled as 

eTmT, efTmfT, efmf. 

The common feature of these topological paramagnets 
is that they all admit gapped surface states with decon- 
fined Z 2 gauge theories, with topological quasi-particles 
labeled as {l,e,m, e}. Notice that we use small letters 
e and m to label anyons on the surface, which are not 
related directly to the bulk E and M particles. For the 
eTmT state, both the e and m particles are Kramers 
doublets with 7^ = — 1. For efmf state, both e and m 
are fermions. The efTmfT state can be viewed as the 
combination of the previous two states, in which both 
the e and m are Kramers fermions. The key property of 
these Z 2 topological order is that they cannot be realized 
in any strictly two dimensional system while preserving 
T. Hence they are called “anomalous”. 

The three topological paramagnets are distinct and 
nontrivial states when existing on their own. But do they 
still give distinct states when combined with a C/(l) spin 
liquid? Or equivalently, is it possible to trivialize the 
corresponding surface topological order in the presence 
of various charged matter fields? 

The stability of topological paramagnets in the pres¬ 
ence of charged matter fields has been studied for some 
cases. It is known that all the topological paramag¬ 
nets are stable if the charged matter field is Kramers 
(T^ = —1), or if the 17(1) charge is T-odd (magnetic- 
like). In the following we show that the eTmT phase 
becomes trivial when the electric charge is either (a) a 
non-Kramers fermion or (b) a Kramers boson. These 
corresponds to EfMi,, EbTMb and E^rMf in Table I. 

The argument is simple: in the presence of an elec¬ 
trically charged particle that is either a non-Kramers 
fermion or a Kramers boson, one can combine that par¬ 
ticle with the e and e particle in the eTmT topological 
order. This is essentially a relabeling of the same phase. 
The resulting topological order is eCmT, which means 
the e particle has charge-1 but is non-Kramers, while the 
m particle is Kramers but charge-neutral. This topolog¬ 
ical order turns out to be realizable even in strictly two 


® The meaning of ‘combining’ is to start from two subsystems, one 
realizing a t/(l) spin liquid, the other realizing an SPT state, and 
couple the two systems weakly. 
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dimensional systems. Hence it is anomaly-free. One way 
to realize this state is to start from the eCTeCT state, 
which is anomaly-free since the m particle is trivial, and 
then put the e particle into a 2D topological insulating 
band. The resulting state is well known [49] to be the 
eCniT. 

The efmf state is nontrivial even in the presence of 
charged particles. The easiest way to see this is to notice 
that the T-broken surface will have nontrivial thermal 
hall conductance, which cannot be canceled by charge 
matter field without introducing another hall conduc¬ 
tance for the gauge field. 

Therefore the topological paramagnets give 
rise to four distinct states when combined with 
EjtMi, and {EfxMf) 0 , and only two 
distinct states when combined with the other three 
phases in Table I. The total number of phases is thus 
4x4-k3x2 = 22. 

IX. PARTON CONSTRUCTIONS 

We now use the insights obtained in previous sec¬ 
tions to describe parton constructions of these C/(l) 
spin liquids. First let us recall that EhMj, is ac¬ 
cessed through the gauge mean field theory, the 
EtTMt through the Schwinger boson representation, and 
EfxMb through the Abrikoson fermion representation. 
Further {EfTMf)g is accessed in the Abrikosov fermion 
representation by putting the fermionic spinons in a 
topological band insulator. 

To obtain “natural” parton constructions for the other 
phases let us consider a spin-1/2 magnet on some lattice 
and use the Abrikosov fermion representation: 

Sr- = (8) 

Here fra is a fermion of spin a =t, i at sites r of the lat¬ 
tice. As is well known this representation introduces an 
SU{2) gauge redundancy and correspondingly the phys¬ 
ical Hilbert space of the microscopic spin system is ob¬ 
tained by imposing a constraint [2]. 

It is convenient for some of our discussion to work 
with Majorana fermions r]ara (ci = 1,2) rather than the 
complex fermions fra- We therefore define 

fra — 2 iVlra f^2ra) (9) 

Let p^,py,p^ be Pauli matrices acting in 771,772 space. It 
is easy to check that the physical spin operators can be 
written as 

S. = ^r7/(pV,n^pV^)77, (10) 

The SU{2) gauge redundancy of the fermion represen¬ 
tation is generated by the operators 

T = ^??‘l77r 


with I = p^, py, <j'^p^). From these we can construct 
SU{2) gauge transformations Or which rotate the Ma¬ 
jorana fermions: 

rjr —> OrPr ( 12 ) 

We first review how to obtain a t/(l) spin liquid 
through this fermionic parton construction before ex¬ 
plaining how to implement time reversal. We consider a 
mean field ansatz described by a Hamiltonian quadratic 
in the fermion operators: 

^mean ~ ^ ^ Pr^rr'Vr' ( 18 ) 

rr' 

with hrr' is a pure imaginary 4x4 matrix. Further 
we must have = —hr'r- Under the SU{2) gauge 
transformation, hrr' gets replaced by Olhrr'Or' ■ 

The unbroken gauge structure is determined by con¬ 
sidering the “Wilson loop” starting from some base point 
r 

Wr[C]=\{hr,r^^, (14) 

The right side is an ordered product over the h-matrices 
connecting the points r^, r^+i that define the closed curve 
C. To get a U (1) spin liquid, all the WrC (for different C 
and bases r) must be invariant under a t/(l) subgroup of 
the full SU{1) gauge group and only this U{1) subgroup. 
In this case there will be a C/(l) subgroup of the gauge 
transformation 

M{ct>) = Y[Mr{ct>) (15) 

r 

(with each Mr{4’) describing an SO{2) gauge rotation by 
angle </>) which leaves the mean field invariant: 

Mr{(j))*hrr'Mr'{(f>) = hrr' (16) 

Let Ur be the infinitesimal hermitian generators of AAfor 
each r (rir will be a purely imaginary, antisymmetric 4x4 
matrix) . Upon including fluctuations Ur will correspond 
to the generators of 1/(1) gauge transformations of the 
spin liquid. 

Now let us consider implementation of physical global 
symmetries. 

In line with the rest of the paper we consider sys¬ 
tems where time reversal is a good global symmetry. We 
make no assumptions about spin rotation symmetry. To 
discuss time reversal properties it is important to distin¬ 
guish two distinct microscopic situations. The physical 
Hilbert space at each site consists of two states - these 
may correspond either to a Kramers doublet or to a non- 
Kramers doublet. Note that this distinction should not 
be confused with the time reversal properties (Kramers 
or not) of the emergent E particle excitations. When the 
physical on-site Hilbert space corresponds to a Kramers 
doublet the spin operators transform under time reversal 


( 11 ) 


as 


S'/ ^ -s/, s+ ^ -s, 


(17) 
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In contrast if the physical on-site Hilbert space corre¬ 
sponds to a non-Kramers doublet we take the spin oper¬ 
ators to transform under time reversal as 

-S^, S+ S- (18) 

For clarity we will focus henceforth on Kramers spin 
systems (Eq. (17)). It is straightforward to extend the 
discussion to non-Kraners spins. Let us now implement 
time reversal on the fermion operators. We may gener¬ 
ally write 

r]r -> frjr (19) 

where T is a 4 x 4 real matrix. Clearly we also have the 
freedom to gauge transform the fermions as part of the 
symmetry implementation, i. e we can multiply T by any 
gauge rotation Or- 

For Kramers spins satisfying Eq. (17) we can take 
fr —>■ fr- This is equivalent to 

f = ia^p^rir (20) 

If a mean field ansatz is time reversal invariant then 
we must have 

f\-hrr')f = Orhrr'Ol, ( 21 ) 

for some gauge transformation Or- The (—) sign in the 
left side is because time reversal is anti unitary and hrr' 
is pure imaginary. Thus we can dehne a “physical” time 
reversal transformation Tr (for any given mean field) 
through 

Tr = fOr ( 22 ) 

under which the mean field is manifestly time reversal 
invariant: 

T^hrr'Tr' = -hrr' (23) 

Now let us consider the algebra of the C/(l) gauge gen¬ 
erators Ur and the physical time reversal transformation. 
The gauge charge W at site r is 

Nr = plnrTjr (24) 

Under time reversal we have 

T-^NrT-^ = r]lT*{-nr)TrVr (25) 

The (—) sign in the right side is because is pure imag¬ 
inary. Generally we have 

T~^nrTr = irir (26) 

The (—) sign describes the group 1/(1) x and the (-I-) 
sign the group C/(l) x Zj (note that if Ur —>■ —Ur, then 
the gauge charge W —> Nr so that the gauge charge is 
even under time reversal). In the former case we should 
take the fermions to be the E particle of the gauge the¬ 
ory, and in the latter we should take them to be the M 


particle. Thus the same parton framework naturally de¬ 
scribes both classes of phases where U is a fermion or 
where M is a fermion. ® 

Note that Tf = TOrTOr- Even though = —1 we 
do not a priori know anything about T^. However we 
know that if is even under Tr, then Tf = El (while 
if Ur is odd we can define a modified time reversal TrMr 
and {TrMrY' can have any value). Thus by choosing the 
mean Held ansatz (which enables us to define Tr and Mr) 
we can access phases where E is either a Kramers singlet 
or Kramers doublet fermion. 

We now use this framework to construct examples of 
the EbMf, {EfTMf)g, and E^tM/ phases. For concrete¬ 
ness we specialize to the three dimensional cubic lattice. 
Consider a mean field ansatz where there is a nearest 
neighbor hopping t and a singlet pairing A on the body 
diagonal. This corresponds to 

hrr' = trr'P^ + Arr'tJ^P^ (27) 

It is easy to check that the fermion spectrum is gapped. 
Further the non-trivial Wilson loops are proportional to 
a^p^ so that this is a gapped 17(1) spin liquid. Corre¬ 
spondingly we have Mr{(j)) = where = -fl 

on the A-sublattice and —1 on the B-sublattice. The 
physical time reversal operator can be simply taken to 
be T = icjyp^ - Thus the generator p^ of Mr is odd un¬ 
der T, and the fermions should be identified with with 
the M particle. Further it is also readily checked that 
the band structure is not topological. We thus have a 
realization of the EgMj phase. 

Next let us modify the t and A to get a topological 
band structure. Precisely such a modification was dis¬ 
cussed in Ref. [50], and requires changing the sign of 
the t and A on some of the bonds. This yields a n = 2 
topological insulator with two surface Dirac cones. The 
17(1) gauge structure and time reversal properties are 
not affected by this modification. We thus end up with 
the Ef,TMf phase. 

Finally, to construct the {EfTMf)^ phase we use 
a different implementation of time reversal. We take 
rjr —> ^rVr Corresponding to Tr = tr- We take a band 
structure in which the two 77 ^ fermions (which make up 
the complex fermion /^) have different dispersion than 
the f fermions. Specifically we choose the f band struc¬ 
ture described in Ref. [50] for the n = 1 topological 
insulator while for the we choose a trivial dispersion^. 
It is easy to check that this mean field ansatz describes a 
17(1) spin liquid, and further that the Ur are odd under 
time reversal. Thus the fermions must be identified with 
M. Further as we have a net n = 1 topological band 
structure we get the {EfTMf)g phase. 


® It is easy to show that if Ur is even under time reversal at one 
site it must be even at all other sites that are connected to it 
and vice versa. 

^ This band structure is invariant under the chosen time reversal 
operation 
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We emphasize that this construction is totally different 
from the “standard” one where the fermions are treated 
as Kramers doublet E particles with topological band 
structure. Nevertheless we get the same phase. 

A. Kramers spin on non-bipartite lattice 

1. Mf-no-go 


rj fermions, the only possibility is = ±1. Clearly 
llb'[C'] must be real to preserve time-reversal. But as we 
discussed above in Sec. IX A 1, on a non-bipartite lattice 
there must exist imaginary Wilson loops. Therefore such 
a construction is impossible. 

Again we emphasize that this does not rule out the 
EfMb state, but does make it less natural in non- 
bipartite Kramers spin systems. 


The previous examples of My-type parton construc¬ 
tion, in which the fermions are monopole-like with their 
[/(I) gauge charge odd under T, were constructed on a 
bipartite lattice (cubic lattice). We now show that for 
Kramers spins (Eq. (17)) such Mf-type construction is 
impossible on a non-bipartite lattice. First we note that 
in Eq. (26) we must choose the -I- sign in this case. Sec¬ 
ond we notice that T = is itself a gauge rotation. 

Thus the physical time reversal matrix = TOr is also 
just an SU{2) gauge rotation. 

Thus we may write = e*®’’”'', where is a hermitian 
generator. In general Tr is a combination of the three 
generators in Eq. (11), satisfying = 1 and r* = —r^. 

With the -I- sign in Eq. (26), we have \Tr,nr] = 0. 
This is possible only if = Err- For any Wilson line 
with base r (see Eq. (14)) satisfies 

Ml{(j))Wr[C]Mr{(j)) = Wr[C] (28) 

As Tr corresponds to a special value of (f) we also have 

TlWr[C]Tr = Wr[C] (29) 

However using Trhrr'Tr' = —hrr' we can also conclude 
that for a loop of length L 

T*Wr[C]Tr = {-l)^Wr[C] (30) 

We thus conclude that L must be even which is possible 
only if the lattice is bipartite. A related argument using 
the trace of the Wilson loop to diagnose time-reversal 
breaking was presented in Ref. [51]. 

This shows that for Kramers spins on non-bipartite 
lattices (such as the pyrochlore), fermionic monopole 
does not arise within the particular (although most com¬ 
mon) type of parton construction from Eq. (8). We 
should emphasize that this does not rule out the pos¬ 
sibility of having such phases in this situation, since one 
can imagine having more complicated types of parton 
construction. However this does suggest that states with 
fermionic monopoles are less natural in these systems. 


2. Ef-no-go 

Following the same logic, we now show that A/-type 
parton construction, in which the fermions are electron¬ 
like but non-Kramers, is also impossible for Kramers 
spins on non-bipartite lattices. To have = 1 on the 


X. LOOP WAVEFUNCTIONS 


It is interesting to understand the differences between 
these different states in terms of their ground state wave 
functions. To that end it is useful to think of the 1/(1) 
spin liquid in terms of fluctuating loop configurations. 
As matter fields (the E and M particle) are gapped, the 
low energy physics is described by Maxwell electrody¬ 
namics. The emergent electric and magnetic fields are 
divergence-free and hence the corresponding field lines 
form closed loops. To describe the wave function we can 
choose either the electric picture or the magnetic pic¬ 
ture (these are different bases for the low energy Hilbert 
space). 

In the specific context of quantum spin ice the mag¬ 
netic flux loops are very easy to picture. Indeed the spin 
ice manifold is parametrized in terms of closed loop con¬ 
figurations formed by the directions of the microscopic 
spins on the pyrochlore lattice. Quantum effects intro¬ 
duce fluctuations of these magnetic loops and, in the spin 
liquid, lead to tensionless fluctuating loops in the ground 
state. The simplest possibility is that the wavefunction 
of the fluctuating magnetic loops is positive definite: 


\^)=Y,'^oiC)\C), 

c 


MC) 


— f d^xd^x'oc 

e 


B(x).B(x') 


(31) 


where a is a positive consitant and B is the magnetic 
field corresponding to the magnetic loop configuration 
C. The positive weight d>o is needed to satisfy Maxwell’s 
equation. Such a “featureless” wavefunction would de¬ 
scribe the EfyMf, phase, as studied in many previous 
works. The M particle is the open end of the magnetic 
loops, and the E particle is a point defect with an addi¬ 
tional phase factor in the wavefunction: 


|4-(E)) = ^e*^^^^/2^o(C)|C), (32) 

c 

where I lables each loop in the configuration C, and 12/ 
is the solid angle spanned by the loop / with respect to 
the E particle. 

To describe the other six phases in Table I, more sub¬ 
tle structures are needed in the loop wavefunction. For 
the EbMf phase, the monopoles-end points of the mag¬ 
netic loops-need to become fermions. This can be done 
by thickening the magnetic loops into “ribbons”, and 
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assigning a phase (—1) to the wavefunction whenever a 
ribbon self-links. More precisely, the wavefunction can 
be written as 

|vI/)=^(-l)^?vi/o(C)|C), (33) 

c 


where is the self-linking number, defined to be the 
linking number of the two boundary loops of each mag¬ 
netic ribbon. An argument in Ref. [52] shows that due 
to this extra phase the open end points of such loops 
have fermi statistics. 

To understand some of the other phases described in 
Table I in terms of fluctuating loops it is more convenient 
to use instead the ‘electric’ picture: the ground state 
is then a superposition of oriented loops (which repre¬ 
sent the electric field lines) with weights derived from 


the Maxwell action d>o ~ e . In the 

EbMb phase these electric loops are featureless, and the 
superposition has positive definite weights for all loop 
configurations just as in Eq. (31). 

In the phase, we can think of the electric field 

lines as ‘stuffed’ with Id Haldane/AKLT chains. One 
way to do it is to consider an additional spin-like order 
parameter n in the disordered paramagnetic phases, and 
assign a Wess-Zumino phase factor in the wavefunction: 


l^f) =^e*’^["('^^'d'o(C)|C;n). (34) 

c 


The easiest way to picture the Wess-Zumino term is to 
view lT[n(C)]/27r as the total skyrmion number of n on 
the membranes whose boundaries are the electric loops 
C. For a closed loop this internal structure has no se¬ 
rious effect. However if we produce an electric charge 
we expose an open end of the electric field line. The 
Kramers doublet known to be present at the open end 
of the Haldane chain then leads to the Kramers degener¬ 
acy of the electric charge. Notice that this possibility is 
meaningful because the electric loop configurations are 
time-reversal invariant. In contrast, the magnetic loops 
cannot be stuffed with Haldane chains, in line with the 
discussions in the rest of the paper. 

If instead the E particle is a fermion (as in EfMb 
or EfTMb) then the electric field is best thought of as 
a thin ribbon (i.e a line with some small but non-zero 
thickness). Again we assign a phase (—1) to an electric 
field loop which has an odd self-linking number, which 
converts the E particle into a fermion. For the EfxMb 
phase, in addition, these electric loops must be stuffed 
with Haldane chains. 

In all the examples above, at least one of the E and 
M particles is trivial (bosonic and non-Kramers). This 
makes it simpler to describe the wavefunction by consid¬ 
ering the loops with nontrivial open ends. For example, 
if the E (M) particle is trivial, we can write the wave- 
function as fluctuating M (E) loops and demand that 
they have the right structure to produce nontrivial end¬ 
points, which are the M (E) particles. 


The remaining two phases in Table. I {EbTMf and 
{EfTMf)e), however, cannot be easily understood using 
the above line of thinking because both E and M parti¬ 
cles are nontrivial. The loop wavefunctions for these two 
phases should capture not only the quantum numbers of 
the end points, but also that of the dual particles. Sim¬ 
ilar issue arises if we want to understand the previous 
phases in the dual loop picture. For example, can we 
have a magnetic loop wavefunction for EfrMb state? 

This issue is actually closely related to the surface 
states of the phases: if a C/(I) spin liquid phase nec¬ 
essarily has a nontrivial surface state, one should be 
able to infer it from the bulk wavefunction. The two 
phases EbxMf and {Efx^f)e both have nontrivial sur¬ 
face states as long as time-reversal is kept. For the 
other phases like EjxMb, the surface has to be nontriv¬ 
ial as long as time-reversal is kept and the M particle 
is not condensed on the boundary. Since the M particle 
is naturally not condensed in the fluctuating magnetic 
loop picture, the wavefunction of magnetic loops should 
contain the information of the nontrivial boundary the¬ 
ory. Similar logic also applies to the other phases except 
EbMb- 

Notice that when the wavefunction is written solely 
in terms of closed loops, the matter fields are naturally 
gapped, even on the boundary. Therefore to have loop 
wavefunctions for the above nontrivial cases, we need 
the loop structures to be able to produce boundary the¬ 
ories with gapped matter fields. Since the gapped mat¬ 
ter fields on a nontrivial boundary is necessarily frac¬ 
tionalized, this suggests the bulk wavefunction should 
be described in terms of fractional loops, instead of the 
“physical” loops such as the 27r magnetic loop. 

A similar problem was tackled in the context of 
symmetry-protected topological (SPT) phases using 
what is known as Walker-Wang const ruction [23, 45, 53- 
55]. The essential idea is that when the boundary is 
a gapped topological order, one can have a loop wave- 
function for the bulk, for which the weights are knot 
invariants of the loop configurations generated by the 
boundary topological field theory. With some modifica¬ 
tion, this idea can be used to generate loop wavefunc¬ 
tions of C/(I) spin liquids (other than the simple EbMb) 
in Table. I. In Sec. XA we discuss a relatively simple 
yet interesting loop wavefunction of the EbxMf phase 
as an illustrating example. With different time-reversal 
implementation the same wavefunction can also describe 
two other phases, namely the EbMf in the electric loop 
picture, and EfMb in the magnetic loop picture, which 
we describe in Sec. XB. In Appendix E we discuss 
a slightly different wavefunction that can describe the 
EbxMb and EfxMb phases. The topological Mott insu¬ 
lator {EfTMf)g can also be described through gauging 
its (non-Abelian) Walker-Wang wavefunction[45], but it 
will be quite complicated and not very illuminating, so 
we will omit the discussion in this paper. 
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FIG. 6: Electric-like loop wavefunction of the Eh^Mf 
C/(l) spin liquid: (a) The amplitude changes sign 
whenever a blue loop links with a red one. (b) A 
doubled blue line can be converted to a doubled red 
line and form a doubled two-segment loop, (c) The red 
and blue loops are switched under time-reversal E- 


A. A loop wavefunction for the E^rMf phase 

For the phase, it turns out to be slightly easier 

to describe the wavefunction in terms of electric loops. 
The wavefunction is written in terms of two species of 
oriented loops, labeled as “red” (r) and “blue” (5): 

I'f') - (35) 

where is the mutual linking number between red 

and blue loops. Two additional features are present in 
the wavefunction. First, a doubled blue line can be con¬ 
verted to a doubled red line and form a doubled two- 
segment loop, even though single lines cannot be con¬ 
verted to each other. Second, the red and blue loops 
get switched under time-reversal action E, which is al¬ 
lowed since these are electric-like loops. These features 
are illustrated schematically in Fig. 6. 

To see that the wavefunction described in Fig. 6 in¬ 
deed describes the E^rMf phase, we need to examine 
the excitation spectrum. First of all, it is useful to ex¬ 
amine the bound state of one blue loop and one red loop, 
with opposite directions. We call this the e loop. No¬ 
tice that due to the condition in Fig. 6b, the e loop is 
undirected. We also note that the e loop has a linking 
sign with both blue and red loops. This makes the end 
points of an individual open blue/red line confined, in 
the sense that they cost an energy proportional to the 
length of the open line. The reason is that with an open 
blue/red line, a small e loop surrounding the interior 
part of the line locally behaves as though it is “linked” 
with the line. However, it cannot have a linking sign in 
the wavefunction. This is because one can continuously 
move the small e loop away from the open line till it 
looks “unlinked”. Therefore the local Hamiltonian near 
the interior of the open line cannot be minimized. Thus 


the energy penalty will be proportional to the length of 
the line. The same physics also appears in Walker-Wang 
models. (See also Ref. [56] for a simple and concrete 
model illustrating this.) 

On the other hand, a doubled red loop (or equivalently, 
a doubled blue loop) can be opened with finite energy 
cost, since it does not have linking sign with anything. 
Therefore we interpret the end points of doubled red 
loops as the deconfined E particles. The single red/blue 
loops are “half” electric loops. 

The open end of a bare e line is conhned because it 
has linking signs with the red/blue loops. However, a 
monopole-like defect can be bound at the end of an e 
line to avoid the sign ambiguity and make it deconfined. 
More precisely, we can have a phase factor in the wave- 
function of an open e line where I de¬ 

notes all the red and blue loops, and is the solid an¬ 
gle spanned by a red/blue loop with respect to the end 
points a, b. This phase factor serves as a smooth interpo¬ 
lator between the linking phase away from the line (-1-1) 
and the linking phase near the interior of the line (—1). 
Since the single red/blue loops are interpreted as half 
electric loops, such a phase factor corresponds precisely 
to a magnetic monopole with unit magnetic charge. 

Therefore the magnetic monopole is bounded to the 
end point of an e line. But notice the e loop is really a 
ribbon, with a red and a blue loop being the edges of 
the ribbon. Therefore it has a self-linking sign from the 
red/blue mutual linking sign, which makes the end point 
a fermion. We have thus obtained fermion statistics of 
the magnetic monopole! 

We now discuss time-reversal action on the E particle, 
which is the end point of a doubled blue line. For this 
purpose it is convenient to view the doubled blue line as 
the combination of a blue, a red and an e line. Under 
time-reversal E the e line is invariant, but the blue and 
red lines are exchanged (see Fig. 6c). Since the time- 
reversed wavefunction away from the charged particle is 
locally indistinguishable from the original wavefunction, 
the E action can be effectively “localized” around the 
E particle, by exchanging the two end points of the red 
and blue lines. Therefore performing E twice amounts to 
twisting the blud-red ribbon, which gives a (—1) phase 
in the wavefunction. This implies that the E particle 
has = — 1 and is a Kramers doublet. 

Similar to the Walker-Wang construction, the above 
discussion is closely related to a possible surface state 
of the EbrMf phase that is gapped but breaks no sym¬ 
metry. This is exactly the M-wall state discussed in 
Sec. VHC. 


B. Alternative loop wavefunctions for EbMf and 
EfMb 

The same wavefunction described in Fig. 6 can also 
describe two other phases, with different ways of imple¬ 
menting time-reversal symmetry. 
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FIG. 7: The same loop wavefunction as in Fig. 6 , but 
with different time-reversal actions, (a) Time reversal 
keeps both the color and the direction of each loop, 
which describes Ei,Mf in the electric picture, (b) Time 
reversal keeps the color, but invert the direction of each 
loop, which describes EfMh in the magnetic picture. 

If time-reversal T keeps both the color and the direc¬ 
tion of each loop (Fig. 7a), then the loops are electric¬ 
like. The argument in Sec. XA for the fermionic 
monopole still applies in this case. But the electric 
charge - the end point of a doubled red (or blue) line 
- now transforms trivially under time-reversal. Hence 
we obtain the E^Mf state in the electric loop picture. 

Now if time-reversal keeps only the color, but in¬ 
vert the direction of each loop (Fig. 7b), the loops are 
magnetic-like. The argument in Sec. X A for the fermion 
statistics of the e particle still applies, but now this par¬ 
ticle should be interpreted as the electric particle E. We 
therefore obtain the EfMb state in the magnetic loop 
picture. 


XI. QUANTUM PHASE TRANSITION 
BETWEEN TWO DISTINCT U{1) SPIN LIQUIDS 

In this section, we consider quantum phase transitions 
from one C/(I) spin liquid to another in Table I. In gen¬ 
eral, one may expect most of the transitions to be first 
order or to go through an intermediate phase, due to 
lack of any obvious order parameter. Thus any contin¬ 
uous transition between two such phases would be quite 
exotic. The understanding of these spin liquids in terms 
of gauged SPT states shed some new light on this sub¬ 
ject. If we can understand possible continuous transi¬ 
tions between different SPT states, we can then couple 
these critical theories to 17(1) gauge fields and under¬ 
stand continuous transitions between 17(1) spin liquids. 
The common feature of such phase transitions is that 
one particle (say E) is unchanged across the transition, 
while the dual particle (say M) drastically changes its 
properties such as statistics, value, and dyon charge. 

Continuous transitions between SPT phases in free 
fermions are well understood [57]. Interestingly, by gaug¬ 
ing such transitions we can already obtain many novel 
transitions between various 17(1) spin liquids. All such 
critical theories are described by massless Dirac fermions 
in (3-1-1) dimensions coupled with a 17(1) gauge field, but 
the effect on the dual particles are very different. We are 
thus able to provide remarkably simple descriptions of 


some highly non-trivial continuous phase transitions be¬ 
tween distinct 17(1) spin liquids. For instance we provide 
a theory for a continuous phase transition between the 
topological Mott insulator and the E^Mf phase. In the 
conventional picture of the topological Mott insulator as 
a spinon topological insulator, such a transition seems 
to require a change of statistics of the electric charge! 
Such a ‘statistics’ changing quantum phase transition is 
however very simply understood within the dual picture 
of the topological Mott insulator (as a monopole topo¬ 
logical insulator) developed in this paper. 


A. “Statistics-changing” quantum criticality: 

Phase transitions of the topological Mott insulator 

1. Warm-up: To EfTMf, 

It is useful to first understand the phase transition 
from the topological Mott insulator - the {EfTMf)g 
phase - to the EfT^b phase. This will set the stage for 
the more surprising (from the conventional viewpoint) 
phase transitions studied below. Since {EfT^f )0 can 
be viewed as a gauged version of a topological insulator 
formed by the Efp particles, we can access its transi¬ 
tion into the EfTMg phase by gauging the topological- 
to-trivial insulator transition. The critical theory is sim¬ 
ply a massless QED with one flavor: 

a^] = CMaxweiiia^i], (36) 

where 70 = U, 7 ^ = aiT 2 ^ V’ = *' 0 ^ 7 o, and time rever¬ 
sal acts as TipT~^ = i(T2'4> (here Ui and are Pauli 
matrices). One can easily check that the term is 
the only T-symmetric mass term. Under a proper UV 
background, to > 0 gives a trivial insulator, which after 
gauging becomes the EfxMg phase, and to < 0 gives 
the topological insulator, which after gauging becomes 
the {EfTMf)g phase, to = 0 is thus the critical point. 

However, to make the transition really continuous, we 
should also forbid other T-invariant terms that close the 
fermion gap before the mass to becomes zero. There are 
two such terms: and . The former is sim¬ 

ply the chemical potential term, and the latter can be 
viewed as a chemical potential alternating in sign for the 
two Weyl fermions. The chemical potential is forbidden 
since the total gauge charge should be zero (assuming the 
fermion gap does not close elsewhere). But to forbid the 
term, more symmetry is required in the theory. 
The simplest possibility is to demand inversion symme¬ 
try X'!/;(x)I“^ = Ti'(/i(—x), as was shown in Ref. [57] 

In general, various anisotropy terms are also allowed. 
These include the spatial anisotropy, the velocity differ¬ 
ence between different bands, and the difference between 
the fermion velocity and the speed of the emergent pho¬ 
ton. We will not go into those details. 
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2. Topological Mott insulator to EbMf 

As discussed in previous sections, the {EfTMf)g phase 
can also be viewed as a “dual” topological insulator 
of the fermionic monopoles. This makes it possible to 
access its transition to the EgMf phase, in which the 
fermionic monopoles form a trivial insulator. The crit¬ 
ical theory has the same Lagrangian as Eq. (36), but 
with a different implementation of time-reversal symme¬ 
try: = iT 2 'lp^ ■ 

Again we need to forbid T-invariant terms that 
can make the fermions gapless. There is only one 
such term, namely the alternating chemical potential 
fjf. We can again forbid it by having inversion 
symmetry x). 

3. Topological Mott insulator to EbxMf 

The previous transition into EgMf phase can be 
viewed as a transition of the fermion monopoles from 
an n = 1 band to an n = 0 band. Now if we consider a 
transition from an n = 1 band to an n = 2 band, this 
gives a transition from the topological Mott insulator to 
the EbrMf phase. The gapless part of the critical the¬ 
ory is the same as the previous transition described in 
Sec. XIA 2. The only difference is that here we have an 
extra gapped Dirac fermion with negative mass, which 
gives the n = 1 band in the UV background. This UV 
background is important in determining the nature of 
the phases away from the critical point. But it will not 
affect the physics right at the critical point such as the 
scaling relations. One can see this by integrating them 
out, which gives a 6 l-term in the U{1) gauge field. But 
since the dual particle (charge) is gapped at the critical 
point, the d-term is a purely surface term and will not 
affect the bulk dynamics. 


B. Kramers-changing quantum criticality: Phase 
transition between EbrMf and EbMf spin liquids 

The EbrMf and EbMf spin liquids differ by whether 
or not the E particle is a Kramers doublet. We now 
discuss the criticality associated with a continuous tran¬ 
sition where this Kramers-ness changes. The EgpMf 
phase can be viewed as fermionic monopoles in an n = 2 
band. We have discussed the transition from the EgpMf 
phase to the topological Mott insulator in Sec. XIA 3. 
We now discuss the transition to the EbMf phase, which 
is a trivial insulator (n = 0 ) of the fermion monopoles. 

The critical theory has two mass Dirac fermions cou¬ 
pled with a U{1) gauge field: 

Maxwell i 

s^l,2 

( 37 ) 


where 70 = n, 7 i = aiT 2 , tps = iipljo, and time re¬ 
versal acts as = iT 24 ’l (here ai and Ti are 

Pauli matrices). Again we need inversion symmetry 
= Ti^/’s(—x) to forbid the alternating chemi¬ 
cal potential term Moreover, we need ini'ijjaips 

to be the unique mass term. This requires some kind of 
rotation symmetry between the two flavors s = 1,2. Mi¬ 
croscopically this could be achieved with certain lattice 
symmetries. We will not go into such details. 

XII. ROLE OF TIME REVERSAL 

We should emphasize that the seven phases are dis¬ 
tinct if and only if time-reversal symmetry is kept. In 
the absence of time-reversal, of course has no mean¬ 
ing, and the 0 -angle can be continuously tuned to any 
value. Even different statistics of particles do not dis¬ 
tinguish phases: if the magnetic particle is a fermion, 
the electric particle is necessarily a boson, and one can 
change the 0 -angle by 27r to shift the m particle to a 
boson[2I]. Similar argument applies if the electric parti¬ 
cle is a fermion. 

Thus in the absence of any symmetries we have exactly 
one C/(I) liquid phase. 

XIII. RELATIONSHIP WITH OTHER WORKS 

We have focused on time reversal symmetric U{1) 
quantum spin liquids. We now place our results in the 
broader context of research on symmetry implementa¬ 
tion in other long range entangled phases. The best un¬ 
derstood long range entangled phases have a gap to all 
bulk excitations, and are characterized by the concept of 
topological order. In the last few years the realization 
of symmetry in such topologically ordered phases has 
gotten a great deal of attention. Such states have been 
dubbed ‘Symmetry Enriched Topological’ (SET) phases. 
As is well known, in an SET phase the topologically non¬ 
trivial quasiparticles may carry fractional quantum num¬ 
bers. This means that the action of symmetry on these 
quasiparticles is non-trivial (technically the symmetry is 
realized projectively rather than linearly). The projec¬ 
tive realization is allowed since these quasiparticles are 
non-local objects. Symmetry operations may also have 
more dramatic effects: they may even interchange two 
different topological sectors. 

In d = 2 space dimensions there is significant progress 
in classifying and understanding these SET phases. For 
some representative papers see Refs. [59-69] and [19, 20, 
55]. For the three dimensional systems of interest in 
this paper, some progress in understanding gapped time 
reversal symmetric Z 2 quantum spin liquids has been 
reported in Ref. [70]. The 17(1) spin liquids discussed in 
this paper are gapless but nevertheless we have shown 
how we can classify and understand the realization of 
time reversal symmetry. 
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We have already discussed previous model construc¬ 
tions of some of these phases. Microscopic models for 
EhMh are common[3-9]. What abont the other phases? 
Ref. [37] constructed a rotor model in which one of the 
emergent particles is a fermion. This can be viewed as 
either a construction of EfMb or of Ei,Mf depending on 
how time reversal is implemented on the microscopic ro¬ 
tor degrees of freedom. If we take it to be a model for 
EbMf it should be possible to modify the Mf hopping to 
give it topological band structure. This will enable writ¬ 
ing down microscopic models for {EfTMf )Q and E^rMf. 
We will however not pursue this here. 


XIV. DISCUSSION: MODELS, MATERIALS, 
AND EXPERIMENTS 

We now consider the lessons learnt from our results for 
current and future possible experimental realizations of 
U{1) spin liquids. We discuss two separate issues. First, 
for a given system, if a 17(1) spin liquid arises, which of 
the seven families of phases in Table I is realized? This 
is particularly relevant to the quantum spin ice materials 
such as the pyrochlore Yb 2 Ti 20 -j. Existing theoretical 
work[ll] assumed that the simplest phase in the E^Mf, 
family is the prime candidate in such systems. However 
this is justified only deep in the spin ice limit[5], and 
materials such as Yb 2 Ti 20 'j are quite far away from this 
limit (see review in Appendix A). So it is important to 
ask which of the seven phases discussed in this paper is 
more likely to arise in such systems. We address this 
issue in Sec. XIV A. 

The next important issue is to identify distinguish¬ 
ing features of these different phases that can be probed 
in experiments. We partially address this issue in 
Sec. XIV B. 


A. Pyrochlore spin ice 

Here we will focus on Kramers spin systems since 
they tend to be more robust against disorder. It is of 
course very hard to decide energetically which phase is 
more favorable, due to the complexity of the underly¬ 
ing Hamiltonian. But at least we can ask the following 
question: which of the seven phases in this paper have a 
natural mean-field description on the pyrochlore lattice? 
Clearly the gauge mean field theory (gMFT) proposed 
in Ref. [II] is a natural mean field theory for the phase 
EbMb- So what about the other six phases? 

As discussed in Sec. IX A, for Kramers spins on a non- 
bipartite lattice such as the pyrochlore, it is quite unnat¬ 
ural - at the level of parton mean field theory - to have 
fermionic monopoles or non-Kramers fermionic electric 
charges. This is already enough to render unlikely the 
EfMb, EbMf and EbrMf phases. 

The phase EbxMb is also unnatural at the mean field 
level on a non-bipartite lattice. To construct it using 


mean field, we need to use the Schwinger boson decom¬ 
position Si_i = The mean field theory should 

have one boson per site on average, and the bosons need 
to be gapped. On a non-bipartite lattice this is impos¬ 
sible at the mean field (quadratic) level, and the bosons 
will always tend to either condense or pair condense, 
which breaks the C/(I) gauge symmetry. 

We are thus left with only two phases: the EfTMb 
and {EfTMf )g. Both can be described at the mean field 
level through the Abrikosov fermion decomposition = 
where time-reversal acts as / —>■ and 

the 17(1) gauge symmetry is the phase rotation on the 
fermions: fa —>■ fa- 

We can try to write down mean field band structures 
of these fa fermions on the pyrochlore lattice. We re¬ 
strict to mean field Hamiltonians that has only nearest- 
neighbor terms, which is reasonable since the spin ex¬ 
change is very short-ranged in quantum spin ice mate¬ 
rials. We further restrict to mean field Hamiltonians 
that are manifestly invariant under the full lattice sym¬ 
metry. This is less justified, since the fa generally are 
allowed to transform projectively under the symmetry 
group. Nevertheless we make this assumption in order 
to make progress, while keeping this caveat in mind. 

With these restrictions, the mean field Hamiltonian 
has only two parameters: the trivial hopping term ti 
and the spin-orbit coupled hopping term t 2 and takes 
the form 

HmF = - flitl + it2drr' ' Cr) fr> (38) 

(rr') 

Here (!„' is a unit vector parallel to the opposite bond 
of the tetrahedron containing rr' (see Ref. [71] for de¬ 
tails). The resulting band structure has been studied 
elsewhere, for example in Ref. [71]. The amusing fact 
about this Hamiltonian is that as long as it is gapped, 
the fermions will always form a topological insulator! 
We therefore reach the conclusion that, beside the EbMb 
state described by gMFT, the topological Mott insula¬ 
tor {EfTMf )0 is the only state that has a simple mean 
field description on the pyrochlore lattice with Kramers 
spins! 

Including fluctuations will lead to a lattice gauge the¬ 
ory with a C/(I) gauge field described by 

the Hamiltonian 

H = -Y,Prtrr'E''^^’fr’ (39) 

{rr') 

(with the hopping matrix trr' = + R2drr' • cr) supple¬ 

mented with the Gauss law constraint 

Y,E„' = flfr-l (40) 

r' 

where E^' is the integer electric field conjugate to o^r'- 
A guess for a spin Hamiltonian which may favor the topo¬ 
logical Mott insulator state is obtained by performing 
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a strong coupling expansion of this lattice gauge the¬ 
ory. The resulting Hamiltonian takes the same form 
as Eq. (Al) for symmetry reason (after the appropriate 
standard rotation from the global basis of spin quan¬ 
tization axis to the local basis). The relative magni¬ 
tude of the coupling constants depend on the parameter 
w = 12 /ti as follows: 

Jzz/J = -|- 8w — 1, 

J± / J — 2rc^ -j- ‘2w -\- —, 

J±±/J = — 2w + 1, 

J^±/J =^/2{-2w^ + W + 1), (41) 

where J is an overall constant, and the underlying 
fermion partons are gapped (and topological) when w < 
—2 or re > 0 and ru 7^ 1. A spin Hamiltonian was ob¬ 
tained in Ref. 31 starting from a Hubbard model. The 
free fermion part of the Hubbard model in Ref. 31, in 
the limit of strong on-site spin-orbit coupling, corre¬ 
sponds to the same hopping Hamiltonian in Eq. (38), 
with w = t 2 /ti = 0.215 determined by orbital physics. 
The spin Hamiltonian obtained in Ref. 31, after a ba¬ 
sis rotation (see Ref. 10), agrees with Eq. (41) with 
w = 0.215. This simple consideration can provide a use¬ 
ful guide in searching for realizations of the topological 
Mott insulator in the family of rare earth pyrochlores if 
their exchange parameters can be determined by experi¬ 
ment. We should caution however that our arguments in 
this subsection are only suggestive, and a reliable deter¬ 
mination of the phase diagram of spin models for these 
pyrochlores is currently beyond the reach of theoretical 
technology. On the experimental side, for Yb2Ti20T, a 
determination of the exchange parameters was provided 
in Ref. 14. However this has been disputed by newer 
experiments [72] which suggest instead a rather different 
set of parameters. In view of the existing uncertainties in 
both the experiment and the theory we will leave further 
discussion of models and materials for the future. 


B. Experimental signatures 

Here we offer some suggestions on experiments that 
may help distinguish these different C/(l) spin liquids. 

Eirst we ask about distinctions in neutron scattering 
experiments. Spin flip excitations that can scatter neu¬ 
trons are created by local operators. If the only global 
symmetry is time reversal, then neutrons will couple to 
all local operators that are odd under time reversal. Let 
us consider a few important ones. Eirst since the emer¬ 
gent magnetic field is time reversal odd, neutrons can 
couple directly to the fluctuations of the internal mag¬ 
netic field. As discussed in Ref. [11] this enables neutron 
experiments to detect the emergent photon. Second the 
number density of magnetic monopoles, and the current 


of emergent electric charge are also local T-odd opera¬ 
tors. Coupling to these will lead to an increase in the 
scattering cross-section when the energy transfer exceeds 
twice the gap of the M and the E particles. Generally 
there should be two thresholds in the scattering cross- 
section set by the M and E gaps. If the E particle is 
a Kramers doublet (which it is in some of the 7 phases) 
then spin flip excitations can also form out of a pair of E 
particles through a combination like E^aE. In this case 
the E- threshold may be much more sharply defined than 
in the case where A is a Kramers singlet. 

Other useful information can be gleaned by studying 
the effects of an applied magnetic field B. Consider 
the phases where A is a Kramers doublet and has a 
gap smaller than M and other composite excitations. 
For a Kramers doublet E particle, a direct coupling at 
quadratic level is allowed: 

AH r-zT3B,E^ajE, (42) 

where is a tensor consistent with lattice symmetries. 
This implies that with increasing B the gap of the E 
particle will close at some finite value before the M gap 
closes. If E is a boson, it will condense in such a B field 
and the U{1) gauge field will be gapped. Since A is a 
Kramers doublet, such a condensate necessarily breaks 
time-reversal symmetry. Therefore the resulting phase 
is magnetically ordered and can be probed through neu¬ 
tron scattering. If A is a fermion, a fermi surface will 
emerge beyond the critical field, and the system becomes 
a [/(I) spin liquid with spinon fermi surface, which can 
be probed through heat capacity or heat transport mea¬ 
surements (see Refs. [73, 74] for interesting recent heat 
transport measurements on pyrochlore magnets). Sim¬ 
ilar consideration also applies if the M gap is smaller 
than the E gap, and can indicate the statistics of the M 
excitation in that case. Thus the behavior in a magnetic 
field can provide useful information to partially distin¬ 
guish these different U{1) spin liquids. 

The two phases {EfTMf )e and E^rMf necessarily 
have protected surface states. As we described it is likely 
that these surfaces are in gapless phases in which case it 
may be possible to detect the surface excitations. A use¬ 
ful experiment will be to deposit a ferromagnet on the 
surface and measure the resulting thermal Hall effect. 
This kind of experiment might perhaps be interesting to 
explore in Yb 2 Ti 20 i. 

XV. CONCLUSION 

In this paper we have provided a detailed understand¬ 
ing of time reversal symmetric U{1) quantum spin liq¬ 
uids. Our results were summarized in Sec. H. To con¬ 
clude we highlight a few open questions. We have not 
discussed the effects of spatial symmetry at all. This 
will lead to a finer distinction between these spin liquid 
phases, and may impact the discussion of existing exper¬ 
imental candidates. A discussion of the effects of space 
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group symmetry on time reversal symmetric C/(l) quan¬ 
tum spin liquids in a cubic lattice was provided in early 
work by Ref. 7 (see also subsequent related work in Ref. 
75 on a model without time reversal on a pyrochlore lat¬ 
tice). Even with just time reversal it will be useful to 
identify sharp experimental fingerprints to distinguish 
the different phases. It may be interesting for future nu¬ 
merical work to study the loop wave functions described 
in Sec. X and explicitly demonstrate their correctness in 
describing the various spin liquids. 

Since the appearance of this paper on the arxiv, there 
have been a number of related further developments 
which we briefly summarize here. Our own subsequent 
work[47] exploited the bulk duality of the topological 
Mott insulator described in this paper to provide a new 
‘dual Dirac liquid’ description of the surface of spin-orbit 
coupled electronic topological insulators. This bulk du¬ 
ality, and the same dual Dirac liquid was independently 
obtained in Ref. 48 which appeared simultaneously with 
Ref. 47. See also the very recent paper [58] with more 
results on the bulk duality. 

We thank Lucile Savary for useful discussions. This 
work was supported by NSF DMR-1305741. This work 
was also partially supported by a Simons Investigator 
award from the Simons Foundation to Senthil Todadri. 


Appendix A: Model for Yb 2 Ti 207 


Here we briefly review the spin Hamiltonian in 
Ref. [10], obtained by fitting the neutron scattering data 
in Ref. [14] through spin wave theory. The values of 
the parameters in this model have, however, been ques¬ 
tioned in more recent work[72]. The Hamiltonian has 
only nearest-neighber terms and takes the form 

= E(.,) - j± (s+s- + S-S+) 

+ j±± {j.,s+s+^:^s-s-) 

+ J.± [S! + QS-) + (* o j)] ,(A1) 

where S'f are spin coordinates in the local basis of spin 
ice. Cijjlij 9-^6 4x4 matrices acting within each tetra- 
hedra and have the form 


c 


/ 0 -1 

-1 0 

giTr/S g—iir/S Q 

, ^-iTr/3 gi-Tr/S _ 


g-iTr/S 

gi7r/3 


"o' / 


,7 = 


The coupling constants are, in meV, 


-c- 


(A2) 


= 0.17 ±0.04, J± =0.05 ±0.01, 

J±± =0.05 ±0.01, =-0.14 ±0.01. (A3) 


If Jzz donimates over the other coupling constants, 
the system at low energy will be restricted to the spin ice 
manifold. But it is not clear whether the above Hamilto¬ 
nian falls into such a regime, given that the other terms 


seem comparable to Jzz in magnitude. Very recent ex¬ 
perimental work[72] has disputed these parameter values 
- the revised values are substantially different and place 
the system even farther away from the regime where the 
restriction to the spin ice manifold is legitimate. 


Appendix B: Time-reversal action on electric charge 

In general, the fundamental electric particle E can be 
multi-component with an internal index i. All the Ei 
particles carry the same gauge charge, which means any 
object of the form EjOijEj must be gauge neutral and 
hence corresponding to a local operator. 

In general, time-reversal could act on E particles as 

TE.T-^ =T,,Ej, (HI) 

where T is a matrix. This implies that 

T^E,T-^ = {T*T\jE,, (B2) 

where the anti-unitarity of T was used. However, on 
any local operator should be trivially identity. Therefore 
we should have 

(T*T)^0{T*T) =0, (B3) 

for any matrix O. This can be true only if r*r = 
where / is the identity matrix. Combining with the com¬ 
plex conjugate relation TT* = we conclude that 

gi0 ^ Therefore T^ET-"^ = ±E. 

Notice that the above derivation assumes that all 
the local objects have = 1 on them, which is true 
for a spin system. However, if charge-neutral Kramers 
fermion is present in the microscopic system (for exam¬ 
ple in a superconductor), then we do have local objects 
with 7"^ = — 1. Following the above logic, it will be pos¬ 
sible to have = ±i for fractionalized objects such as 
E particles. 


Appendix C: The same U{1) spin liquid from 
gauging two distinct insulators 

We start from fermions with C/(l) x E symmetry, 
which corresponds to the fermionic monopoles discussed 
in the main text. It is known that for free fermions the 
band structures of such kind of fermions are classified by 
an integer n. We now show that the 17(1) spin liquids 
obtained by gauging the fermionic monopoles with band 
topology n is the same as that from gauging the fermions 
with band topology —n. The reason is very simple: at 
the level of band structures, the difference between a ±n 
band and a —n band is the way time-reversal T is imple¬ 
mented. This can be most easily seen through the surface 
Dirac cones, where the effective Hamiltonian looks the 
same for bands at ±n: 

n 

H = '^1p\ {Px(^x ± PyCrz)i’i, (Cl) 

7=1 
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but time-reversal acts differently for ±n: = 

dsziayip'^. 

Before gauging, the different time-reversal action 
makes it impossible to continuously tune one state into 
the other. But after gauging, the difference becomes sim¬ 
ply a gauge transform U = where Q is the charge. 
Therefore the ±n bands become gauge equivalent, and 
give rise to the same quantum spin liquid. 

It is also known[25, 26] that a band indexed by n la¬ 
bels the same interacting bulk state as a band indexed 
by n -|- 8. It is also known that bands indexed by n and 
n -|- 4 differ from a bosonic symmetry-protected topolog¬ 
ical state called eTmT topological paramagnet[19, 20], 
which is protected only by time-reversal. Together with 
the previous identification of n and —n bands, we con¬ 
clude that states with band index n = ± I (mod 8) and 
n = ±3 (mods) corresponds to two distinct phases. The 
bulk excitations of these two phases are identical and 
correspond to the {EfTMf)g U{1) spin liquid (topologi¬ 
cal paramagnet). The two phases have different surface 
states, and one can obtain one state from the other by 
combining with a eTmT topological paramagnet. 

As discussed in Sec. VIII, the eTmT topological para¬ 
magnet becomes trivial when combined with the EhxMf 
state. Therefore all the states with band index n = 
2(mod4) corresponds to a unique E^rMf state. 


Appendix D: Surface states of various 17(1) spin 
liquids 

Here we discuss the surface states of the first five 
phases in Table I. All these phases have a particle {E 
or M) that is a trivial boson. Therefore the “wall” cor¬ 
responds to this particle is trivial. However, they can 
still have a nontrivial wall corresponding to the nontriv¬ 
ial quasiparticle. For example, the EgMf phase could 
have an M-wall, with a Z2 topological order {1, e, m, e} 
on the wall, where the e and m particles carry electric 
charge qe = 1/2. Following the logic in Ref. [21], this 
wall will convert the fermionic monopole into a bosonic 
one upon tunneling. The other U{1) spin liquids have 
similar nontrivial walls, with a Z2 topological order and 
proper symmetry assignments on the e and m particles. 

The physics of these “walls” will be important when 
we put two different spin liquids side by side. For exam¬ 
ple, if we put an EgMg next to an EgMf, can we have 
coherent tunneling between the quasiparticles across the 
interface? Naively this is impossible since the monopole 
is fermionic in one region but bosonic in the other, and 
the system will just be two 17(1) spin liquids essentially 
decoupled from one another. However, we can put the 


A-wall described above on the interface, and the M 
monopoles can now tunnel through the interface, with 
an e particle left on the wall. 



FIG. 8: Loop wavefunctions of the EgTMg, EfTMi, and 
other phases: the amplitude changes sign whenever a 
directed loop links with an undirected one. 


Appendix E: Loop wavefunctions of various U{1) 
spin liquids 

Here we discuss loop wavefunctions of various other 
C7(I) spin liquids in the same spirit of Sec. XA. One 
can think of this approach as starting from a Walker- 
Wang wavefunction[23, 45, 53-55] and “gauging” the 
C7(I) symmetry. Many of these loop wavefunctions can 
be written in a form slightly different from that in 
Sec. X A. The wavefunctions have a directed “half” loop 
and an undirected loop condensing simultaneously, with 
a (—1) phase in the wavefunction whenever a “half” loop 
and an undirected loop mutually link (see Fig. 8). Fol¬ 
lowing similar analysis in Sec. XA, the resulting phase 
is a 17(1) spin liquid, with two fundamental particles E 
and M. One of them is bound with the open end of 
an undirected loop, and the other is the open end of a 
doubled half loop. 

If the directed loops reverse their directions under 
time-reversal, they are magnetic loops; otherwise they 
are electric loops. It is now easy to see how to assign 
various quantum numbers to the dual particles in this 
description. For example, in the magnetic loop picture, 
we can put a Haldane chain in the undirected loop to 
make the E particle Kramers, in which case we obtain 
the EbTMg phase in the magnetic loop picture. If we 
also make the undirected loop a ribbon with a (—1) self¬ 
linking sign, we convert the E particle to a fermion and 
obtain the EfxMb phase in the magnetic loop picture. 

There are some more complicated cases, including the 
Eb^Mf phase in the magnetic loop picture, and the 
{EfTMf )e in both pictures. The surface topological or¬ 
der of these phases are well studied in the literature, from 
which one can derive the corresponding Walker-Wang 
wavefunctions and their gauged version[23, 25, 26, 43- 
46]. However, the results are somewhat complicated (es¬ 
pecially for {EfTMf)g which is non-Abelian) and not 
particularly illuminating, so we will obmit the discus¬ 
sion here. 
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